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of  the  imager- tracker  system. 
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SUMMARY 

This  technical  report  describes  the  results  of  a  two-year 
investigation  into  the  Advanced  Imaoinq  Tracker  (AIT),  a  novel 
infrared  tracking  concept  which  attempts  to  provide  low  resolution 
imagery  at  frame  rates  comparable  to  contemporary  FLIR's  with  simul¬ 
taneous  high-bandwidth  target  centroid  estimation.  The  initial 
motivation  was  to  find  useful  alternative  passive  tracking  systems 
for  the  DARPA  Talon  Gold  and  LODE  programs. 

The  program  was  divided  into  two  phases.  Phase  I,  a  concept 
feasibility  study,  was  completed  in  December  1980.  The  results  of 
this  phase  were  combined  in  an  It  erim  Report.  The  Phase  II  pro¬ 
gram  effort  centered  on  fabricating  a  workinn  AIT  portable  bread¬ 
board,  containing  all  optics,  electronics,  processing  hardware  and 
software  needed  for  conducting  an  AIT  experimental  test  program. 

The  breadboard  was  constructed  and  AIT  imaqing  and  tracking 
demonstrated.  During  Phase  II  new  classes  of  imaging  and  tracking 
algorithms  were  defined,  analyzed,  characterized  by  computer 
simulation  and  implemented  in  hardware.  A  special  purpose  hiyh-speed 
digital  signal  processing  system,  the  Programmable  Microcoded 
Processor  (PMP),  was  designed,  fabricated,  tested  and  coded.  The 
optical  system  is  capable  of  visible  and  IR  operation,  and  includes 
a  dynamic  target  simulator  for  system  testing.  An  extensive  soft¬ 
ware  base  was  written  to  support  algorithm  development,  system 
simulation,  experimental  data  analysis,  and  interactive  control. 

A  color  display  system  was  integrated  with  the  breadboard  system  for 
pseudo-color  image  output  and  diagnostic  display.  The  system  is 
ready  for  use  in  an  experimental  program  of  AIT  imaging  and  tracking 
performance. 
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1 . 0  INTRODUCTION 
1  •  1  AIT  Concept 

The  Advanced  Imaging  Tracker  (AIT)  program  undertook  to  develop 
a  novel  tracker-imager  concept.  The  scheme  employs  a  quad-cell 
detector  and  an  image  scanning  device  to  extract  from  a  received 
optical  image  its  centroid  and  a  low-resolution  reconstruction  of 
the  image.  The  system  was  originally  intended  as  an  advance  on 
scanning-r'LIR  type  imaging  trackers,  since  it  could  provide  a  much 
higher  centroid  tracking  rate  (4-8  kHz,  as  opposed  to  30  H2  for 
video-rate  FLIR's)  while  retaining  the  FLIR  image  rate  (30  Hz). 

The  AIT  concept  is  described  in  somewhat  more  detail  in  Figure 
1.1.  Passive  infrared  radiation  (dashed  line)  from  a  distant  target 
is  received  by  a  telescope  and  transmitted  to  the  AIT  through  some 
optical  relay  system  including,  in  this  case,  a  set  of  fast  tilt 
mirrors  for  fine  track  on  the  target.  In  this  example,  the  telescope 
is  part  of  an  HEL  beam  director  system,  so  that  a  Shared  Aperture 
Device  is  included  in  the  optical  path. 

The  optical  layout  of  the  AIT  itself  is  shown  figuratively  at 
bottom  left  in  the  diagram.  Here  the  major  components  of  the  AIT 
are  shown.  A  carefully  prepared  optical  plane  wave  (solid),  whose 
amplitude  may  be  controlled  by  a  modulator,  is  combined  with  the 
incoming  IR  from  the  telescope,  and  reflect  d  into  the  remaining 
tracker  optics.  The  plane  wave  is  used  as  a  highly  accurate  tilt 
reference,  which  may  be  used  for  differential  measurement  of  the 
image  centroid  in  the  focal  plane,  removing  measurement  errors 
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FIGURE  1.1 


induced  by  motion  of  the  AIT  optical  components,  electronics 
offsets,  etc. 

The  combined  beam  is  passed  through  a  hioh  speed  nutation 
device,  depicted  here  as  a  tip-tilt  mirror.  This  device  deflects 
the  beam  in  a  complex  pattern  in  the  focal  plane,  where  it  is 
detected  by  a  simple  quadrant  detector.  The  result  is  a  set  of 
four  output  wavefroms,  one  from  each  detectoi ,  which  are  passed 
to  a  pair  of  signal  processors. 

The  first  of  these  is  based  on  the  centroid  extraction  tech¬ 
nique  (the  Synchronous  Centroid  Sampling,  or  SCS,  algorithm)  which 
was  developed  by  AOA  for  the  1^  Sensor  wavefront  sensor.  The 
details  of  this  process  will  be  discussed  in  Section  3,  but  its 
purpose  is  to  determine  the  IR  centroid  of  the  image,  and  hence 
angular  position  of  the  target  object. 

The  second  processor  is  used  to  reconstruct  the  target  image 
from  the  detector  wavefroms.  The  image  information,  obtained  at 
relatively  slow  speed,  can  then  be  combined  with  the  high  speed 
centroid  data  to  determine  the  appropriate  aim  point  for  the  HEL 
beam  on  the  target  object.  The  computed  aim  point  results  are  then 
applied  to  the  tip-tilt  mirrors  to  maintain  system  track. 

The  AIT  process  therefore  consists  of  three  basic  steps: 

l)  Nutate  the  image  on  a  quadrant  detector  t  <  produce 
the  output  waveforms  which  carry  the  centroid  and  image  information. 
?.  I  Process  the  waveforms  for  centroid  data. 

3)  Process  the  waveforms  for  image  reconstruction. 
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1 . 2  Aim-Point  Maintenance 

The  need  for  both  image  and  centroid  data  is  outlined  in 
Figure  1.2.  Whereas  the  centroid  of  the  image  may  be  obtained 
at  high  speed  (at  least  using  the  AIT  technique)  it  does  not  provide 
sufficient  information  to  stabilize  an  HEL  beam  on  a  relatively 
small  region  of  an  extended  target.  This  is  because  the  target 
will  in  general  be  thermally  dynamic;  that  is,  its  temperature 
profile  will  change  in  time  due  both  to  its  intrinsic  characteristics 
and  its  interaction  with  the  impinging  HEL  beam.  In  the  example  of 
Figure  1.2,  the  target  is  an  atmospheric  missile.  Due  to  air 
frictional  heating  and  the  different  thrust  phases  of  its  flight, 
the  temperature  distribution  on  the  target  will  change  significantly 
as  the  leading  edges  heat  up,  the  engine  area  warms,  and  the 
exhaust  plume  vanishes  at  burnout. 

Besides  these  intrinsic  effects,  the  target  will  develop  a 
hot  spot  at  the  point  of  initial  contact  with  the  HEL  beam.  For 
a  tracker  system  in  which  the  HEL  aim  point  is  not.  coincident  with  the 
centroid  measuranerit,  the  contact  point  would  be  offset  from  the  centroid  of 
the  IR  emission  of  the  target  alone.  As  the  hot  spot  develops, 
however,  the  centroid  will  be  shifted  towards  the  hot  spot,  which 
will  in  turn  change  the  dim  point.  This  interaction  will  cause 
the  hot  spot  to  drift  on  the  object,  greatly  reducing  the  effec¬ 
tiveness  of  the  HEL. 

If  information  from  the  full  image  is  available  in  addition  to 
its  cont'Oid,  the  HE  l  diinpoint  may  be  calculated  as  an  offset  to 
the  IR  centroid.  This  offset  can  than  be  updated  as  the  target 
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FIGURE  1.2 
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temperature  profile  changes,  stabilizing  the  HEL  aimpolnt.  Since 
the  temperature  profile  is  governed  by  thermal  time  constants,  it 
vaires  slowly,  on  the  order  of  tens  of  milliseconds,  compared  with 
the  sub-millisecond  jitter  rates  of  the  centroid  itself  (due  to 
effects  such  as  atmospheric  disturbance  and  telescope  bearing 
jitter). 

The  AH  is  therefore  intended  to  yield  the  combination  of 
fast  centroid  tracking  and  slower  image  reconstruction  (at  a  ratio 
of  about  100:1)  required  for  aim-point  maintenance. 
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1 . 3  A IT  Program  Goal s 


The  AIT  program  was  designed  to  analyze  the  concept  and  identify 
the  relevant  theoretical  and  technological  issues,  characterize  the 
problems  and  develop  solutions,  and  finally  implement  the  results 
in  a  working  breadboard  imager-tracker. 

The  problems  of  imaging  and  tracking  were  addressed  separately. 
The  tracking  problem  was  treated  as  an  extension  of  the  I  Sensor 
measurement  ter1  ique,  and  the  basic  system  constraints  imposed 
by  the  need  tu  track  were  determined.  The  imaging  problem  was  then 
analyzed  subject  to  the  tracking  constraints,  and  a  formal  proce¬ 
dure  for  obtaining  image  reconstruction  developed.  The  tracking 
and  imaging  algorithms  were  then  evaluated  by  extensive  computer 
simulation  to  verify  the  analytical  results.  A  processing  archi¬ 
tecture  capable  of  imp! ementi nu  the  algorithms  was  adopted,  and  a 
complete  processor  designed  and  fabricated,  and  coded  to  execute 
both  imaging  and  tracking  processes.  A  breadboard  version  of  the 
optical  sensor  head  itself  was  then  constructed,  and  integrated 
with  a  specially  built  target  simulator  to  test  the  performance  of 
a  complete  system.  Finally,  operation  of  all  components  was 
demonstrated,  and  image  reconstruction  experimentally  achieved. 

In  the  following  sections  each  of  the  major  areas  of  program 
effort  a^e  described  and  results  presented.  Th  .  material  is 
organized  topically  rather  than  chronologically,  since  many  of  the 
efforts  were  parallel;  the  time  sequence  implied  by  the  section 
numbers  is  however  roughly  correct. 
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The  AIT  program  was  divided  into  two  parts:  a  Phase  I  feasibility 
study  which  extended  from  June  3980  to  December  1980,  and  a  Phase  II 
breadboard  fabrication  and  test  program  conducted  from  January 
1981  to  January  1982.  The  results  of  Phase  I  have  been  repeated 
and  extended,  so  that  this  report  documents  the  output  of  boch 
phases. 
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2 . 0  THEORETICAL  FORMULATION:  IMAGING 
2 • 1  The  I3  Sensor 

The  AJT  tracker- imager  is  an  extension  of  the  I"*  Sensor ^ 
centroid  tracker,  and  it  is  useful  to  first  examine  the  simpler 

3 

I  Sensor  for  insight  into  the  operation  of  the  AIT, 

3 

Figure  2.1  is  a  scnematic  view  of  the  I  Sensor.  Designed 
as  an  optical  wavefront  sensor,  it  is,  in  its  simplest  form,  an 
array  of  tilt  sensors,  each  measuring  the  mean  t , 1 t  in  a  sub¬ 
aperture  of  the  instrument's  full  input  aperture.  Tne  resulting 
2-0  array  of  wavefront  slope  measurements  can  then  be  used  to 
obtain  ari  estimate  of  the  original  wavefront  shape. 

In  operation,  the  input  wavefront  is  combined  with  a  plane 
wave  which  acts  as  a  tilt  reference.  The  reference  is  switchable, 
and  is  used  only  a  small  fraction  of  the  time  to  calibrate  the 

sensor  for  its  own  optical  and  electronic  drifts.  The  wavefronts 

g 

are  then  relayed  through  a  nutation  device,  in  the  case  of  the  I‘ 
sensor,  which  imposes  a  simple  circular  scan  on  the  beams.  The 
full  aperture  is  then  divided  into  subapertures  by  an  array  of 
sampling  lenslets.  At  the  focal  plane  of  each  lenslet  is  a 
quadrant  detector. 

The  detector  focal  plane  is  shown  enlarged  at  the  bottom  of 
the  figure.  As  the  subaperture  focal  spot  moves  in  the  circular 
pattern  induced  by  the  nutation  scan,  the  detector  quadrant  output 
vary.  The  intensity  in  ^ach  quadrant  increases  as  the  spot  moves 
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into  that  quadrant,  at  the  same  time  the  intensity  in  the  adjacent 
quadrant  decreases  as  the  spot  leaves.  When  appropriately 
demodulated,  as  described  in  Section  3,  the  centroid  of  the  unnutated 
focal  spot  (or  the  centroid  of  the  circular  orbit)  is  obtained  in 
each  subaperture. 

Taken  individually,  each  subaperture  tilt  sensor  is  a  complete 
centroid  tracker.  Under  the  Advanced  Imaging  Adaptive  Optics  (AIAO) 
program,  it  was  shown  that  the  variance  a  on  the  angular  measurement 
is  given  by 

o 2  =  J-p  (^)2  (Radians)2  (2.1) 

where  SNR  =  signal -to-noise  ratio  of  the  full  four-quadrant 
detector  area  in  the  operating  bandwidth 
X  -  wavelength  of  incident  radiation 

D  =  the  subaperture  diameter 

or,  for  photon-noise  limited  operation 

a2  =  f  (^)2  (radians)2  (2.2) 

where  N  =  number  of  photons  detected  during  observation  period 
This  is  within  a  factor  of  two  of  the  theoretical  performance 

(21  'j 

limit  for  quad-cell  type  detection,  v  '  so  that  the  I  sensor  is 
a  near-optimal  centroid  tracker. 
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2.2  Nutation  of  Extended  Objects 

By  examination  of  the  detector  waveforms  for  an  extended 
object,  it  can  be  seen  that  more  information  than  just  centroid 
data  has  been  encoded  by  the  nutated  quad-cell  process.  Figure  2.2 
shows  the  quad  cell  outputs  for  a  horizontal  bar  target  having  a 
5:1  aspect  ratio.  The  waveforms  were  experimentally  measured  on  a 
breadboard  sensor  set-up.  The  oblong  nature  of  the  bar  is  apparent 
by  the  slope  differences  in  the  waveforms  for,  say,  the  Q1-Q2 
transition  and  the  Q2-Q3  transition.  As  the  long  axis  of  the  bar 
crosses  between  quadrants  the  transition  is  gradual,  while  the 
transition  is  abrupt  when  the  narrow  axis  moves  across  the 
quadrant  boundary. 

The  nature  of  the  transformation  between  image  and  nutation 
may  be  quantified  by  considering  the  object  to  be  composed  of 
pixel  elements  of  varying  intensity  in  the  image  plane.  An  illus¬ 
tration  of  the  results  of  such  a  decomposition  is  snown  in  Figures 

2.3  and  2.4.  In  Figure  2.3,  three  different  pixels  from  a  field  of 
8x8  are  shown  illuminated.  Because  of  their  differing  displace¬ 
ments  from  the  quad  cell  origin,  the  arcs  described  by  each  will 
result  in  different  residence  times  for  each  pixel  in  the  four 
quadrants.  Simulated  quadrant  outputs  for  each  pixel  are  shown 

in  Figure  2.4.  From  these  considerations  one  is  invited  to  make  a 
one-to-one  correspondence  between  each  pixel  and  its  characteristic 
set  of  waveforms.  In  the  next  subsection,  the  time  sampled 
versions  of  these  waveforms  will  be  called  pixel  vectors. 


...  .  .  |,  ■ 
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NUTATED  QUAD  CELL.  OUTPUT  FOR  EXTENDED  OBJECT 
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The  process  of  nutation  and  quad  cell  detection  can  therefore 
be  viewed  as  a  linear  transformation  between  the  pixel  represen¬ 
tation  of  the  image  and  a  pixel  vector  representation.  This  concept 
is  illustrated  in  Figure  2.5.  In  order  that  the  pixel  vector 
representation  be  complete,  certain  sampling  constraints  must  be 
observed.  As  will  be  seen  in  subsection  2.3,  this  consideration 
leads  to  the  use  of  more  complex  nutation  patterns  than  the  circle 
(which  likewise  complicates  the  tracking  process,  as  seen  in 
Section  3).  The  image  reconstruction  problem  is  then  seen  as  ana¬ 
logous  to  using  an  inverse  transformation  of  the  nutation  process 
to  obtain  the  original  image.  In  the  next  section  a  procedure  for 
finding  this  transformation  is  developed. 


ui . . 
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2 • 3  The  AIT  Image  Reconstruction  Algorithm 

The  AIT  imaging  algorithm  reconstructs  the  input  image  by 

solution  of  the  matrix  equation  which  relates  the  detector  output 

waveforms  to  the  image  intensity  distribution  for  a  given  nutation 

pattern.  In  order  tc  define  this  equation,  we  consider  an  N  x  N 

element  square  pixel  field  which  will  contain  the  image.  If  the 
t  h 

amplitude  of  the  i  pixel  is  given  by  the  scalar  A.; ,  then  the 
entire  image  can  be  defined  by  a  vector  containing  all  pixel 
ampl i tudes. 

t  =  A, 


A  2  (2.3) 

N 

Consider  an  image  in  which  only  the  ith  pixel  is  illuminated 
with  unit  amplitude,  i.e., 

Ai  -  1 

Ak=i  -  0  (2.4) 

This  will  result  in  a  specific  detector  output  waveform  when 

nutated  on  the  quad  cell.  If  the  four  detector  outputs  are  time 
') 

sampled  with  4N^  samples  taken  as  a  full  nutation  pattern  is 
scanned,  the  1 6N ^  detector  output  samples  corresponding  to  unit 
illumination  of  the  i ^  pixel  define  the  vector,  called  a  pixel 


vector: 
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P.  =  P 


11 


QM 


P44N2 


1-4  Quadrant  number 


M  -•  1  -  4N  Time  sample  number 


(2.5) 


There  is  a  unique  pixel  vector  £  associated  with  each  pixel 
in  the  image  field.  The  detecto'  output  vector  15  is  defined  simi¬ 
larly  to  ?  but  is  the  response  to  a  general  image  input  The 
vector  15  can  be  described  as  a  linear  combination  of  the  unit  pixel 
vectors  f*.  and  this  relationship  can  be  formalized  by  arranging  the 
pixel  vectors  to  form  a  forward  matrix  F 


F  =  P, 


P  2 
N 


(2.6) 


1  '  ' 

-► 

The  relationship  between  the  input  image  A  and  detector  outputs 
D  can  be  stated 

F  A  =  D  (2.7) 

Since  [5  is  the  measured  variable  and  A  the  desired  result 
(the  image),  this  equation  must  be  solved  for  A.  This  can  be  done 
formally  by  defining  a  backward  matrix  B  which  is  an  inverse  of  F. 

B  F  =  I  I  is  the  identity  matrix  (2.8) 


That  gives 


A 


B 


(2.9) 


matrix  backward 
uiatri  x 


detector- 

vector 
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.  p  p 

The  matrix  F  is  non-square  (having  dimension  N  x  1 6 N  for  a 
nutation  pattern  with  ideal  sampling),  so  that  B  can  be  any  of 
the  possible  pseudo-inverses  of  F.  The  pseudo- inverse  having 
minimum  variance  for  noisy  input  is  defined  by 


§»i  „'<-fT£>'V  <2-I0> 

Reconstruction  of  an  image  from  quadrant,  detector  outputs 
requires  a  nutation  pattern  sufficiently  complex  to  unambiguously 
sample  all  pixels  in  the  required  reconstruction  field.  That  is, 
the  spatial  Nyquist  criterion  of  at  least  two  samples  per  linear 
dimension  per  pixel  must  be  obtained.  In  terms  of  AIT,  the  image 
must  be  nutated  so  that  the  quad  cell  axes  fall  at  four  distinct 

sample  points  in  each  pixel.  For  perfectly  uniform  sampling, 

?  ? 
this  implies  4N  samples  per  quadrant  for  a  total  of  16N  samples 

needed  for  an  alias-fee  reconstruction,  where  N  is  the  number  of 

pixels  across  one  edge  of  a  square  field. 

Here  the  spatial  Nyquist  criteria  of  sufficient  sampling  is 

specified  by  the  existence  of  the  inverse  (F1  £) 
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2 . 4  Prope rti  e s  o_f_  1 1 he  Recon s t rue tjo n  A Vgojri  t hm 

As  stated  earlier,  the  back  matrix  B_  (equation  2.10)  will  yield 

(3 ) 

the  minimum  norm  estimate  of  the  image  from  the  input  data  D, 

Tne  exact  relationship  between  the  variance  on  the  reconstructed 
image  pixels  and  the  system  noise  is  derived  below. 

2.4.1  AIT  Imaging  Noise 

The  AIT  imaging  process  is  a  multiplexed  detector 
approach.  Although  the  nutated  quad  cell  detector  arrangement 
yeilds  a  very  efficient  centroid  tracker,  its  use  as  an  imager 
is  a  compromise  from  a  signal -to-noise  standpoint.  In  the 
following  sections  an  expression  for  the  noise  propagator  from 
detector  output  to  image  output  is  derived. 

AIT  Noise  Pro pa gator 

We  will  use  the  definitions: 

A. j  =  amplitude  of  i^  pixel  in  input  image  (i^  entry 

of  image  vector) 

£\A.  -  noise  in  i  u  pixel  (zero  mean  random  variable) 

B. .  back  matrix  coefficient  for  i  p’xel  and  jtri 

i  J 

detector  sample 

0.  -  detector  sample  (j™  entry  of  detector  vector) 

J 

/\D.  =  detector  sample  noise  (zero  mean  random  variable) 
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We  wish  to  determine  the  variance  squared  in  the  i  pixel  after 
one  frame.  This  will  be  given  by 

<(Ai  +AA.J ) 2>  *  <( l  B.  ■  (Dj+ADj  ) ) 2>  (2.11) 

J 

where  <f  )>  denotes  expected  value,  and  the  sum  is  over  the  j 
detector  values. 


Expanding 
<(«,)2  * 
<<VZ> 


0  since  AA^ 
has  zero  mean 


2(1  B„DjI  B.jiOj)  *  (J  > 


<( 


> 


<(AAj)2> 


o,vr’>  ♦  <?();  BijOji  B(/oj)>  ♦ 

<o;  »,iAD/> 

J 


variance  squared 
in  i th  p j xel 


0 


s i nee  AD  .  has 
J 

zero  mean 


=  <  cr  Bfj6DJ)Z>  (2.12) 

=  <0,  BjJADj)> 
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2.4,2  AIT  image  con d it i onj ng  characteristics 

The  formulation  of  the  reconstruction  algorithm  alleys 
a  certain  amount  of  image  conditioning  to  be  included  in  the  back 
matrix  13,  which  can  have  significant,  engineering  implications. 
Specifics! ly  if  the  forward  matrix  F  is  generated  by  a  particular 
hardware  AIT  system  by  sequentially  presenting  the  systen  with 
point  source  illumination  in  each  of  the  pixels  in  the  input 
representation  field  and  recording  the  detector  outputs,  F  will 
contain  all  the  effects  of  the  system  optical  transfer  function. 

The  back  matrix  B  will  include  the  inverse  of  the  OTF.  As  a  result 
resolution  loss  due  to  optical  aberration,  image  distortion,  and 
other  image  defects  can  be  removed  in  the  reconstruction  process. 
The  extent  to  which  this  can  be  done  depends^ ft  the  characteristics 
of  the  system  OTF;  significant  amounts  of  deblurring  have  been 
achieved  using  similar  techniques.  This  image  conditioning 

property  arises  as  a  desirable  side  benefit  cf  this  reconstruction 
approach. 
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3 . 0  THEORETICAL  FORMULATION:  TRACKING 
3 . 1  Algorithm  Constraints 

The  twin  needs  of  imaging  and  tracking  act  as  competitive 
constraints  on  the  nutation  pattern  used  to  scan  the  object.  The 
spatial  Nyquist  criteria  (Section  2.3)  dictate  the  sampling  density 
needed  to  provide  imagery  at  a  specific  resolution.  The  desire  for 
fast  tracking,  using  an  1°  Sensor-like  tracker  algorithm,  requires 
a  sampling  sequence  which  circles  the  quad  cell  origin  many  times 
for  each  complete  sample  set  in  order  to  get  many  centroid 
estimates  per  image  estimate.  In  addition,  engineering  limitations 
on  the  nutation  device  itself,  which  will  be  discussed  in  Section 
7,  have  led  to  use  of  a  specific  class  of  compound  Lissajous 
figures  for  nutation  patterns.  The  image  field  is  deflected 
according  to  the  following  set  of  parametric  equations: 

Xnut  =  Dmax  [(sin2Trfpt)  (cos2Tr(fn/s)t) j  (3.1) 

Ynut  =  aTlax  ((cos2lTfl,t)  (sin2ir(fn/s)t)j 
Umax  =  pattern  size, 
f^  =  nutation  frequency, 

S  =  number  of  spirals  per  pattern. 

This  pattern,  referred  to  as  the  collapsing  ellipse  pattern, 
can  be  seen  as  composed  of  a  circular  nutation  pattern  of  nutation 
frequency  f  ,  multiplied  ly  an  amplitude  modulation  function  at 
frequency  f  /s.  The  sampling  sequence  therefore  makes  S  circuits 
about  the  quad  cell  center  for  every  complete  pattern.  An  example 


wmmmmmmmm  > 
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of  this  pattern,  which  was  used  as  the  benchmark  throughout 
the  program,  is  shown  in  Figure  3.1.  The  pattern  was  designed 
to  cover  an  8  x  8  pixel  imaging  field.  The  parameter  S  15, 
and  there  are  31  sample  points  taken  every  spiral,  for  a  total  of 
31  x  15  =  465  samples  per  image. 

Given  these  basic  nutation  characteristics,  a  tracking 


algorithm  can  be  defined. 


PBr“’  1  '! " 1 1 1  '"ir*"*  ■»*<-’  iiw^ifjfiin 
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3.2  T fie  SCS  Algor i thm 


The  SCS  algorithm  can  be  defined  by  reference  to  Figure  3.2. 

A  focal  spot  corresponding  to  the  target  image  is  nutated  on  a 
quad  cell.  The  quadrant?  are  labeled  -  Q^.  Four  intervals  in 
the  nutation  cycle  are  also  defined,  as  periods  T^  -  T^ ,  or  nutation 
angular  intervals  0^  -  8^,  which  are  equivalent  since  with  circular 
nutation  the  spot  velocity  is  constant.  Measurements  of  spot 
centroid  displacement  are  made  at  different  times  for  the  two  axes-, 
x-displacement  is  measured  during  periods  T,  and  T_,  while  y  is 

a  0 

measured  during  T?  and  T^.  The  algorithm  defines  the  spot  position, 
and  therefore  input  tilt,  by  a  numerator  proportional  to  displace¬ 
ment  and  brightness,  and  a  denominator  proportional  to  brightness 


only: 


XNUM 

XDEN 


(Q1(T2)_Q2(T2))  "  <Q4(T2)“Q3(T2) }  +  (Q4(T4)_(Q3(T4) }  " 

(Q1(T4)‘Q2(T4))  (3-2) 

(q1(t2)+q2(t2))  -  (q4(t2)+q3(t2) )  +  (Q4(T4)+(Q3(T4))  + 

(Q1(T4)+Q2(T4))  {3-3) 


XNUM 

XDEN 


(3.4) 


The  motivation  for  this  definition  can  be  seen  by  separately 
examining  the  effects  of  each  term.  The  differences  Q, (T0)-Q0(T„ ) 
arid  Q4(T4)-Q3(T  >  correspond  to  measuring  the  intensity  imbalance 
between  right  and  left  half-planes,  during  those  times  when  the 
snot  is  expected  in  those  quadrants,  which  occurs  due  to  spot 
displacement.  The  imbalance  occurs  because  a  shift  of  the  spot 
appears  as  a  shift  in  the  center  of  the  nutation  circle,  as  shown 
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in  Figure  3.3.  The  imbalance  between  the  signals  integrated  from 
each  quadrant  can  be  interpreted  as  proportional  to  the  differences 
in  arc  lengths  bounded  by  the  integration  intervals  and  quad  cell 
boundaries,  as  shown  in  the  figure. 

The  other  numerator  terms  Q4 ( Tg ) “Qa ( )  an^  Q] ( T4 ) -Qg ( ^ 4 )  are 
the  power  differences  seen  when  the  focal  spot  is  not  present  in 
the  respective  quadrants.  For  an  ideal,  sharply  bounded  spot,  the 
second  set  of  terms  should  represent  differences  0 n?y  in  background 
radiation  or  detector  responses,  and  serve  to  cancel  these  contri¬ 
butions  from  the  first  set,  so  that  the  numerator  becomes  sensitive 
to  the  spot  only.  For  spots  whose  radii  exceed  the  size  of  the 
nutation  radius,  or  for  less  well  bounded  spots  (such  as  the  Airy 
pattern  or  Gaussian  distributions  commonly  encountered)  the  spot 
intensity  is  never  wholly  absent  from  a  quadrant  even  though 
nutation  may  have  moved  the  spot  farthest  away  from  that  quadrant 
during  a  cycle;  in  these  cases  spot  power  as  well  as  background 
power  is  cancelled  from  the  numerator,  and  modulation  efficiency  is 
reduced.  The  relationship  between  object  size  and  nutation  radius 
is  therefore  an  important  consideration. 

For  objects  smaller  than  the  nutation  radius,  the  scale  of 
the  displacement  is  given  by  the  nutation  radius,  that  is.  the 
output  of  the  algorithm  gives  displacements  as  fractions  of  the 
nutation  radius.  The  response  of  the  algorithm  as  a  function  of 
nutation  radius  is  summarized  in  Figure  3.4.  At  small  radii,  the 
sensitivity  is  determined  by  the  spot  size,  while  at  large  radii  the 
nutation  radius  sets  the  sensitivity. 
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The  function  of  the  denominator  is  to  normalize  the  output 
to  spot  brightness.  It  can  be  seen  as  a  sum  of  left  and  right 
half-plane  terms,  rather  than  a  difference  as  in  the  numerator. 
Background  cancelling  terms  are  still  present. 

The  algorithm  for  y-axis  information  is  similarly  defined, 
using  the  time  intervals  T^  and  T^.  Note  that  y-axis  information 
is  obtained  90°  out  of  phase  with  x-axis  information  during  the 
cycle. 

It  is  important  to  emphasize  that  the  background-cancelling 
terms  make  the  SCS  algorithm  an  AC  algorithm.  The  position 
information  is  impressed  on  a  carrier  at  the  nutation  rate.  Out 
of  band  noise,  in  particular  1/f  noise,  is  rejected  by  the  same 
cancelling  action  as  removes  background  illumination. 
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3.3  The  Alf  Tracking  Algorithm 

Significant  modifications  must  be  made  to  the  SCS  algorithm 
to  accommodate  the  AIT  nutation  pattern.  Some  performance  loss 
might  be  expected,  but  the  following  features  should  be  retained: 

1)  AC  operation.  The  continued  recognition  that  many 
of  the  AIT  applications  will  be  in  infrared  systems  means  that 

AC  demodulation  must  be  implemented  to  avoid  1/f  noise  and  uniform 
background  illumination  problems. 

2)  Bandwidth.  A  closed  figure  is  described  once  every 
15  times  the  spot  orbits  the  center;  imaging  is  updated  every  time 
a  full  pattern  is  completed.  This  will  be  referred  to  as  a  major 
cycle,  as  shown  in  Figure  3.5.  The  separate  traversals  will  be 
denoted  minor  cycles;  one  example  is  shown  in  Figure  3.6.  In 
order  to  ratain  the  AIT  feature  of  tracking  at  greater  rates  than 
imaging,  tracker  data  should  be  obtained  at  the  minor  cycle  rate, 
which  is  15  times  that  of  the  major  cycle,  or  imaging  rate.  This 
complicates  the  problem  further,  since  each  minor  cycle  is  different 
from  the  others,  requiring  a  d i f ferenf.  algorithm  form  for  each. 

3)  Well-defined  scale  factor.  In  the  SCS  algorithm, 
the  scale  factor  is  defined  by  the  nutation  radius  as  1  c rig  as  the 
object  size  is  smaller  than  that  radius.  In  the  AIT  nutation 
pattern,  it  is  clear  that  there  exists  no  single  nutation  radius, 
so  that  the  scaling  of  the  output  becomes  less  jbvious.  Since 
varying  size  objects  will  be  observed,  the  algorithm  must  be 
adaptable  to  varying  object  extent. 
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Using  these  criteris,  we  now  outline  a  general  prescription 
for  an  AIT  tracker  algorithm,  which  uses  a  form  analogous  to 
eq.  (3.21  -  (3.4): 

S"™<M)  ‘  ^I'V'VV’  ‘  <W  -  W>  + 

N4(M)(tW  -  W>  -  N^ll/W  -  W>  (:‘-5) 

Xden(H)  ’  +  W'  -  N2(M)(W  *  W>  * 

C~~  c— 

N4(M^W  +  "  N.^WW  +  ^3'6^ 


X  (M) 

X(M)  =  V%“(TfT  *  Xoff(M) 


(3.7) 


Here  the  quadrant  outputs  are  referenced  to  a  sample  number 
hi  of  a  minor  cycle  M,  where  for  the  current  AIT  N  <  31,  M  <_  15, 
the  separate  sets  of  samples  are  analogous  to  the  intervals 

Tj  -  in  the  SC5  algorithm,  and  define  which  samples  are  used  to 

form  an  interval.  This  procedure  will  be  discussed  in  detail 
below.  The  term  XQf.^(M)  is  ;n  offset  correction  term  that  accounts 
for  residual  asymmetry  in  the  sample  set,  which  differs  for  diffe,  nt 
minor  cycles.  The  factor  S(M)  is  a  scale  correction  factor,  which 
defines  a  nominal  or  average  nutation  radius  for  the  cyc’e. 

Consider  the  example  of  Figure  3.7.  A  minor  cycle  is  shown 
(actually  tne  most  difficult  of  the  rvcle  shapes  tc  treat),  with 
its  31  sample  segments  marked.  The  two  circles  indicate  two 
example  object  sizes  of  radius  1  pixel  \7  pixel  diameter)  and  2 
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pixels  (4  pixel  diameter).  To  measure  an  X-displacernent.  using  an 

SCS-type  algorithm,  the  object  must  be  entirely  above  or  below  the 

x-axis.  For  '.'-measurements,  the  object  must  be  completely  to  the 

left  or  right  of  the  y-axis.  This  insures  that  the  measurement  is 

independent  of  object  size.  Based  on  this  notion,  ail  samples 

outside  Region  1  can  be  included  for  the  smalt  obejct,  and  those 

outside  Region  2  can  be  used  with  the  large  object. 

A  second  constraint  is  that  a  sampling  interval  must  extend 

at  least  an  object  width  to  either  side  of  the  axis  which  it  crosses 

during  the  measurement  interval.  The  intervals  ultimately  defined 

by  these  constraints  are  shown  in  Figures  3.8  and  3.9  for  the  two 

and  four  pixel  diameter  objects  respectively. 

Note  that  for  the  four  pixel  spot,  x-axis  information  cannot 

be  ot  ined  in  this  minor  cycle,  since  there  are  no  samples  in 

Quadrant  1.  f- information  for  this  large  object  would  have  to  be 

obtained  f ;  om  other  minor  cycles  with  more  amenable  shapes. 

The  response  of  the  algorithm  for  any  chosen  minor  cycle  and 

object  size  can  be  calculated  for  a  point-object,  essentially  by 

measuring  the  x-  or  y-axis  components  of  the  included  arc  lengths 

in  the  measurements.  By  obtaining  values  of  X(M)  for  a  zero  offset 

spot  and  a  unit  offset  spot  (for  a  given  M),  the  constants  S(M) 

and  X  cart  be  obtained  from  equation  (3.7). 

oft 
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FIGURE  3.8.  Usable  arc  lengths  for  2  pixel  diameter  object. 
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3.4  AIT  Tracker  S/N 


The  expression  for  angular  measurement  variance,  eq.  (2.1), 
of  the  I  Sensor  can  be  modified  to  provide  an  estimate  of  the  AIT 
tracker  S/N  performance.  Two  observations  must  be  made.  One,  eq. 
(2.1)  is  true  for  a  system  operating  at  optimal  nutation  radius, 
that  is 


where  is  the  nutation  radius,  X  the  received  wavelength,  and 
D  the  subaperture  diameter.  The  ratio  A/D  is  the  angular  resolution 
of  the  aperture,,  so  it  is  seen  that  best  performance  occurs  with 
the  angular  scan  radius  equal  to  the  pupil  angular  resolution. 

In  general,  for  nutation  radii  larger  than  A/D,  the  variance  is 
given  by: 

°2  =  SNR  (RN)2  (3-9) 


where  is  the  (not  necessarily  optimal)  angular  nutation  radius. 

For  AIT,  R^  can  be  replaced  by  a  mean  or  normalized  nutation  radius 
corresponding  to  an  average  over  the  arc  used  in  each  minor  cycle. 

For  AIT,  R^  will  nearly  always  be  larger  than  optimal,  to  maintain 
object  size  insensitivity. 

In  addition  to  the  variance  increase  encountered  due  to  nutation 
size  increase,  there  is  also  a  duty  cycle  reduction  that  results 
from  use  of  only  part  of  the  total  available  arc  length;  the 
loss  to  performance  is  not  exactly  equal  to  the  fractional  loss 
of  arc  length  since  the  samples  taken  at  different  nutation  angles 


41 


contribute  differently  to  the  measurement.  As  a  worst  case, 
however,  one  can  assume  that  the  variance  is  proportionate  to  the 
fraction  of  a  nutation  cycle  used: 

Om 

6-  •  =  —  3.10) 

pi,J  tt  v  ; 

|L  a.  U 

where  3.  •  refers  to  the  i  axis  (x  or  y)  of  the  jtn  minor  cycle, 
and  0^  is  the  total  nutation  included  by  the  minor  cycle. 

An  upper  bound  on  the  AIT  tracker  variance  can  therefore  be 
given  as 


2  2 
SNR  RN  6i  -i 

1  »o 


(3.11) 
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4.0  SOFTWARE  SIMULATION  AND  PERFORMANCE  CHARACTERIZATION:  IMAGING 

All  of  the  image  reconstruction  algorithms  were  characterized 
using  an  extensive  set  of  simulation  codes  run  on  the  AOA  Data 
General  computer  system.  This  section  discusses  the  structure  and 
function  of  the  major  simulation  packages,  the  procedures  involved 
in  constructing  the  back  matrix,  and  the  results  obtained  using  the 
process  to  reconstruct  images  from  simulated  data  sets. 

4.1 _ Computer  System 

The  simulations  were  run  on  the  Air  Force-owned  Data  General 
Nova  3/12  system  located  at.  AOA.  The  conf-iguration  of  this  micro¬ 
computer  system  is  diagrammed  in  Figure  4.1.  The  computer  is 
equipped  with  64K  words  of  mapped  memory,  hardware  floating  point 
capability  and  a  ULM  I/O  port  for  communication  with  extra  external 
devices  (such  as  the  printer  and  the  PMP)  under  RS- 232  protocol. 

The  system  supports  two  interactive  users  in  foreground/barxground 
mode.  A  floppy  disk  and  a  cartridge  disc  drive  (5  MByte  fixed,  5 
removable)  are  for  mass  storage.  Hard  copy  of  text  is  output 
through  a  DG/Dasher  printing  terminal. 

Two  peripheral  systems  communicate  with  the  Nova  through  its 
backplane,  the  DG/DAC  system  and  the  color  display.  The  DG/DAC 
laboratory  interface  includes  4  A/D  channels,  D/A  channels,  and 
16  bits  each  of  latched  TTL  inputs  and  outputs.  This  system  enables 
computer  access  to  experimental  hardware  for  control  and  data 
col  1 ect i on  purposes . 
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This  interface  was  also  used  to  access  the  AIT  interface 
electronic  system,  which  contains  the  nutation  drive  generator  and 
a  switched  integrator  front,  end  for  data  sampling  of  the  nutation 
waveforms. 

The  Genisco  Programmable  Graphics  Processor  and  associated 
hardware  provide  high-speed,  medium  resolution  co'or  graphics  capa¬ 
bility  for  pseudo-color  image  output.  This  particular  system  was 
chosen  ror  its  high  data  transfer  rate  and  color-fill  speed,  which 
could  accommodate  the  anticipated  maximum  AIT  image  rate. 
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4.2  S ijnu la t i on  Structu re 

To  provide  unambiguous  characterization  of  algorithm  performance 
and  retain  the  operational  flexibility  needed  to  conduct  algorithm 
development,  the  program  structure  shown  in  Figure  4.2  was  adopted. 
Depicted  are  three  system  segments,  two  of  which  are  software 
constructs . 

Complete  simulation  of  a  proposed  AIT  configuration  first  requires 

a  means  for  generation  of  the  nutation  waveforms  associated  with 

the  desired  input  test  pattern.  This  function  is  performed  by 

„  1 

program  5IMDSK.  Originally  written  to  simulate  I  Sensor  outputs 
0 

under  the  Sensor  Study  program,  SIMDSK  was  expanded  to  accommo¬ 
date  the  more  complex  nutation  patterns  used  in  AIT.  SIMDSK  is 
designed  to  parallel  as  closely  as  possible  the  physical  process 
as  done  in  the  experimental  set-up.  The  program  accepts  a  test 
image  defined  as  a  20  x  20  array  of  intensity  values.  This  input 
file  is  then  numerically  superimposed  on  a  40  x  40  element  array 
which  defines  the  quad  cell.  The  quad  cell  routine  can  include 
the  effects  of  nonuniform  responsivity  and  finite  gaps  between 
detector  elements.  The  2-D  integral  of  the  power  in  each  detector 
is  calculated  and  output. 

The  displacement  [’'•tween  the  image  and  quad  cell  coordinate 
systems  can  be  set  by  the  operator  (as  an  input  tip  and  tilt)  and 
by  another  file  containing  the  description  of  the  nutation  pattern. 

The  detector  output  calculation  is  repeated  for  each  of  the  spot 
positions  defined  in  the  nutation  file.  The  result  is  a  vector  of 
detector  values  given  as  a  function  of  nutation,  i.e.  a  quad  cell 


waveform  set. 
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The  simulated  cpjad  cell  outputs  a«~e  stored  in  a  disk  file  for 
later  access.  Generation  of  a  detector  waveform  can  take  up  to  an 
hour,  depending  on  the  resolution  used  in  the  calculation. 

For  reconstruction  of  tfie  nutated  image,  the  program  IMAGE  is 
used.  The  input  to  image  is  the  disk  file  containing  a  nutation 
waveform.  A  reconstruction  algorithm  is  then  called,  which  is  se¬ 
lectable  by  the  user.  Normally  this  would  involve  data  conditioning 
steps  and  then  multiplication  of  the  detector  vector  by  the  back 
matrix.  The  output  image  vector  can  then  be  displayed  by  a  variety 
of  peripheral  devices,  such  as  a  printer  plot,  grey-scale  oscillo¬ 
scope  display,  or  the  Genisco  color  display. 
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The  computational  magnitude  of  the  pseudo-inverse  generation 
step  is  considerable,  and  at  this  time  the  problem  will  be  scaled 
by  introducing  the  specific  parameters  used  in  the  AIT  program,  as 
listed  in  Table  4.1. 


Parameter 

Size 

Quantity 

Dimension 

Pixel  field  size 

8x8 

A 

64  pixels 

Number  of  minor 
cycles  per  pattern 

15 

S 

Number  of  samples 
per  minor  cycle 

31 

M 

Number  of  samples 
per  pattern 

4  x  15  x 

?.l  6 

I860 

Size  of  forward 

64  x  1860 

F 

119,040 

matrix 

coefficient 

Size  of  normal 
matrix 

64  x  64 

(ftf) 

4096  entries 

Size  of  back 
matri  x 

1860  x  64 

B 

119,040 

TABLE  4.1.  PARAMETERS  OF  AIT  RECONSTRUCTION  PROBLEM 


Because  of  the  limited  computation  speed  and  particularly  the 
limited  memory  oize  of  the  Nova  computer,  the  problem  had  to  be 
broken  into  many  sections.  Two  kinds  of  routines  were  designed 
One  type  is  a  set  of  computational  routines  for  performing  matrix 
algebraic  operations  such  as  matrix  multiplication  and  inversion. 
The  other  is  a  set  of  disk  file  handling  programs  for  packing  and 
manipulating  the  large  matrices  and  minimizing  time  consuming 
disk  accesses  by  careful  buffering.  Particularly  important  is  the 
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design  of  the  transpose  routine,  which  can  take  up  to  six  hours 
to  transpose  a  forward  matrix  if  no  storage  symmetries  are  exploited. 

The  computational  steps  used  to  generate  a  back  matrix  8  are 
listed  in  Figure  4.3.  Each  block  shown  contains  the  name  of  the 
program  step,  the  name  of  the  main  sub.outine  called  (in  parenthesis) 
and  the  name  of  the  output  file  (following  the  arrow).  The  output 
file(s)  for  eacn  step  form  inputs  for  the  subsequent  steps. 

The  first  version  of  this  code  required  nearly  twelve  hours 
of  computer  time  to  generate  L  it  once  the  reconstruction  para¬ 
meters  were  fixed  and  the  buffe.  mg  optimized,  the  back  matrix 
generation  time  was  reduced  to  about  an  hour. 

During  development  of  these  codes,  several  important  numerical 
constraints  were  discovered,  pertaining  to  calculational  accuracy 
needed  at  various  stages  in  the  procedu  e: 

1.  The  determinants  cf  the  (FTF)  matrices  are  very  large 
'of  order  106<r),  and  can  achieve  values  during  the  inversion  process 
which  overflow  the  machine  if  attention  is  not  paid  to  this  problem. 

?. .  The  three  computation  intensive  steps,  in  which  F  F, 
(FTF)~'*,  and  (F  r)  *F^  are  formed,  must  be  done  using  double  preci¬ 
sion  arithmetic  for  meaningful  results  to  be  obtained  (on  the  16 
bit  Neva  3). 

3.  The  input  data,  i.e.  the  nutation  pattern  NUTXY  and 
the  forward  matrix  F  need  not  have  double  precision  accuracy  in 
order  to  obtain  a  well-defined  solution.  This  has  important,  impli¬ 
cations  for  the  hardware  impelmentation,  suggesting  that  there  will 
be  no  unanticipated  restrictions  on  the  accuracy  of  the  nutator  or 
data  collection  process. 
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I L  [  PS  ,  (N'JTL  I  PS )  -  NUTXY 
Define  nutation  pattern 


GrMX ,  ( FORW )  - >  F  .MX 

Define  forward  matrix  F 


GFTF,  (MMULT  ,DMULT)-»-  FTF .  MX 
r0ru  FTF 


?DINV,(MINV) - >  FTF.INV.I1X 

Perform  inversion  (_FTF)~* 


|  MFINV 

j  Verify  (FTF)(FTF)4  -  l 

L _ L_ _ -'I _ _ _ 


TRANSP.(MTP)  -*■  FT.TX .TFNV.TX 

Transpose  matrices  to  facilitate 
mu1 ti pi ication 


GBACK.(MHULT)  -  BACK. MX 

Multiply  ( FTF )" 1  F1  -  B 

SBACK 

TBACK 


FIGURE  4.3 
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♦  •  4  linage  Jjuto  ;  t,y 

The  ev.pected  quality  of  the  images  produces  by  this  system  can 
be  deduced  by  examination  of  t>e  identity  matrix  formed  by  the  pro¬ 
duct  B  F  =  1.  We  ask  the  quesf  >n:  how  weii  are  single  pixels 
reproduced  by  the  imaging  process? 

The  first  row  of  the  identity  generated  by  the  minimum  variance 
B  F  product  is  shown  in  figure  4.4.  1 ne  B-matrix  was  obtained  using 

toe  procedure  described  in  the  previous  section,  with  the  F-matrix 
generated  by  the  15/31  collapsing  ellipse  pattern.  "he  row 
contains  the  coupling  coefficients  for  each  of  64  'nput  pixels  into 
output  pixel  number  one.  The  coupling  is  truly  excellent.  Input 
pixel  one  is  transferred  to  output  pixel  one  with  unit  amplitude 
to  the  accuracy  of  the  representation.  Contribution  from  other 
pixels  to  output  pixel  one  are  down  by  at  least  1?  orders  of 
magnitude.  Similar  accuracy  is  obtained  for  the  63  other  pixels. 

As  a  point  of  comparison,  toe  first  row  of  the  B  f  product 
resulting  from  a  buck  matrix  generated  by  an  iterative  technique 
previously  developed  during  the  program  (see  Phase  I  Interim 
Report)  is  shown  in  Figure  4.5.  The  coupling  to  the  correct  pixel 
is  only  about  93%,  while  significant  transfer  occurs  from  other 
input  pixels  (note  especially  pixel  10  -  row  2  ei  !umn  2  -  which 
has  a  coefficient  of  -  126.  This  is  a  very  significant  Dixel- 
pixel  interaction).  These  results  were  still  subjectively  acceptable, 
as  seen  in  the  PHASt  I  interim  report  results. 
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FIGURE  4.5 


4.5  I  val  uati  on  of  I ma cj e  No i  s e 


Equation  2.15  in  Section  2.4.1  defined  the  equation  for  the 
noise  associated  with  the  AIT  reconstruction  process: 


<  (aa.)2> 

1  AIT 


y  B.  .2  (P  )2 

j  ^  ^FRAME  NAI1 


(2.15) 


Once  a  back  matrix  has  been  defined,  the  quantity  T.B.j  can 
be  evaluated.  For  comparison  purposes,  expressions  for  the  pixel 
noise  associated  with  other  imaging  techniques  have  been  derived: 


For  a  staring  array  FlIR 


V  B.  .  =  1  ,  since  only  one  coefficient  B . .  =  1 

4  1J  U 


Ny  =  1  (no  scanning) 


<I(AA  )^>  = _ -1 

^»SA  AT 


^NPIXEL^ 


FRAME 


(4.1) 


P^PIXel  is  noise  power  for  a  single  array  element. 


For  a  single  detector  scanning  FLIR 


Nt  =  NTcn  (number  of  steps  for  a  two-dimensional  scan, 
1  minimum  of  N). 


<  (AAi )  >  $0 


Vn(P 


NPIXEL ' 


FRAME 


(4.2) 
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For  a  linear  scanning  FLIR 


Ny^<.  (number  of  steps  for  a  one- dimensional  scan, 


minimum 


of  /N 


V  B.j  -  1,  since  only  one  non -zero  coefficient 


«m,F>,_s 


N  t  p  ) £ 

WTLS'  NPIXF!./ 


AT 


FRAME 


(4.3) 


In  order  to  make  a  direct  comparison  between  the  effects  of 
the  imaging  techniques  themselves,  comparable  imaging  conditions 
must  be  defined.  First,  the  frame  time  '•>  set  to  one. 

Then  the  relationship  between  the  detector  noise  power  for  AIT  must 
be  correctly  written.  To  cover  a  specific  field  of  view,  each 
pixel  is  mapped  into  a  corresponding  solid  angle  of  target  space. 

The  mapping  of  angular  f.o.v.  in  target  space  onto  the  corresponding 
detector  area  is  illustrated  in  Figure  4.6.  It  is  seen  that  for  a 
given  field  angle  0,  the  detector  size  d  required  for  proper 
coverage  is 


2D0 

NA 


(4.4) 


where  D  is  the  telescope  aperture  diameter  and  NA  is  the  numerical 
aperture  of  the  detector  condenser  optics. 

If  0  is  the  total  imaging  f.o.v.,  then  in  v.he  FLIR  systems, 
each  detector  covers  the  angular  region  corresponding  to  one  pixel. 
For  an  N  x  N  field,  the  pixel  field  angle  is  0/N,  or 
2  DO 


FUR 


(4.5) 
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FIGURE  4.6.  RELATIONSHIP  BETWEEN  ANGULAR  F.O.V.,  AND  DETECTOR  AREA 


56 


For  AIT,  each  detector  quadrant  must  cover  the  entire 
f.o.v.,  so  the  AIT  detector  size  is  given  by  equation  (4.4). 

Since,  for  IR  detectors,  the  detector  noise  power  is  given  by 

9 

P NOISE  =  dlf<AGERB  D*  (4.6) 

where  P  is  the  signal  bandwidth,  and  D*  is  the  normalized  detectivity 
of  the  detector. 

Comb:ning  equations  (4.4)  -  (4.6),  we  obtain 

%AIT  ~  ^  '  PNFl IR  (^etector  or  background  noise  limit)  (4.7 

In  the  photon  noise  limit,  the  noise  is  proportional  to  the 
area  Af  of  the  target,  rather  than  the  area  of  the  detector,  so 

PNAIT  =  At  '  PNFLIR  (Photon  noise  1iinlit)  (4-3) 

The  imaging  signal  to  noise  performance  of  idealized  FLIRs 
versus  several  versions  of  the  AIT  are  collected  in  Table  4.2. 

The  comparison  assumes  an  3  x  8  pixel  field.  The  "15/31  Ellipse" 

AIT  is  the  system  studied  under  this  program,  with  a  mean  noise 
factor  of  1.05  obtained  by  averaging  the  noise  factors  of  all  the 
6’  pixels.  The  "ideal"  AIT  assumes  a  perfect  raster  scan  nutation 
pattern  tc  minimize  noise  factor,  which  is  estimated  at  .25.  This 
kind  of  pattern  would  be  unsuitable  For  tracking. 

Performance  for  photon  noise  limit  can  be  obtained  by  the 
substitution  indicated  in  equation  (4.3). 
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COMPARATIVE  S/N  PERFORMANCE  OF  DIFFERENT  IMAGER  TYPES 

NOISE  POWERS  VARY  AS  DETECTOR  AREA. 

AIT  DETECTORS  MUST  COVER  FULL  F.O.V., 

FLIR  DETECTORS  COVER  ONLY  SINGLE  PIXEL 
F.O.V. 

FOR  N  «  64  PIXELS 

PNAIT  “  ^  PN PIXEL  "  8  PNPIXEL  =  8  PN 

ATFRAME  "  1 

/ \  Bij  VnOISE^ 

J 


IMAGER 

#  OF 

DETECTORS 

'fi* 

nt 

P NOISE 

-  /N  ^  a 

P 

STARING 

ARRAY 

64 

1 

1 

PM 

V 

LINEAR 

SCAN  FLIR 

8 

1 

8 

PK 

2.82  °N 

SINGLE 

EL E MEN f 
FUR 

1 

1 

64 

PN 

3PN 

AIT  (15/31 
ELLIPSE) 

4 

1.05 

1860 

3pN 

353.5  PM 
N 

"IDEAL" 

AIT 

4 

(.25) 

?56 

8Pm 

64  PN 

TABLE  4. 2 


The  reconstruction  matrix  can  be  tested  for  image  formation 
using  the  combination  of  programs  S1MDSK  and  IMAGE.  The  input  to 
the  SIMOSK  program,  which  generates  the  simulated  waveforms  is 
defined  as  a  numerical  intensity  file.  An  example  is  shown  in 
Figure  4.'/,  which  is  the  input  file  for  a  2  x  4  pixel  horizontal 
bar,  located  at  coordinates  (Jv  0,  and  having  unit  intensity. 

The  output  of  SIMDSK  is  a  file  containing  the  nutation  waveforms 
corresponding  to  the  input  image.  Hgure  4.8  shows  a  plot  of 
the  waveforms,  unpacked  to  show  the  outputs  of  quadrants  1-4 
individually. 

These  output  waveforms  can  then  be  used  as  inputs  to  the  IMAGE 
program  whi-.h  multiplies  the  detector  outputs  by  the  reconstruction 
matrix  and  displays  the  results.  The  pseudo-color  output  for  the 
simulated  horizontal  bar  is  shown  in  Figure  4.9. 

The  results  shown  are  essentially  perfect.  This  is  due 
partly  of  course  to  the  high  accuracy  of  the  reconstruction 
algorithm,  but  also  to  the  fact  that,  the  input  object  was  defined 
in  perfect  registration  with  the  pixel  field,  so  that  no  effects 
of  spatial  quantization  of  the  image  are  encountered. 

A  mere  realistic  case  is  the  reconstruction  of  the  same  bar 
rotated  45°,  as  shown  in  Figure  4.10,  Here  die  effects  of  spatial 
quantization  are  obvious,  and  the  results  are  more  representative 
of  image  coding  at  this  coarse  resolution. 
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The  effects,  of  quantization  become  more  apparent  in  reconstruc¬ 
tion  of  continuous- intensity  objects,  figure  4.11  shows  a  small 
Gaussian  profile  spot  reconstruction,  and  Figure  4.12  shows  a 
square  with  Gaussian  smoothed  edges.  As  well  as  the  coarseness 
of  reconstruction,  an  axis  of  diagonal  symmetry  is  apparent  in 
what  are  defined  to  be  circular  and  four-fold  symmetric  objects. 

This  feature  is  an  artifact  of  the  sampling  pattern,  which, 
especially  near  the  center  of  the  field,  has  non-uniform  sampling 
reflecting  the  same  diagonal  symmetry. 


64 


4 . i  Algorithm  Parameter  Sen >si tivi ty 

The  simulation  codes  were  used  to  test  the  sensitivity  of  the 
reconstruction  algorithm  to  small  variations  in  system  parameters 
of  engineering  importance.  Tested  were  sensitivities  to 

1)  Waveform  shift.  Routines  were  written  to  numerically 
snift  the  detector  waveforms  in  time,  this  was  done  to  simulate 
the  effects  of  phase  shifts  in  real  detector  circuits  and  phase 
error  in  synchronism  between  the  nutation  scan  and  sampling 
ci  rciPtry. 

?.)  Waveform  smoothing.  Simulates  the  effects  of 
band  1  ini i ti ng  in  the  readout  electronics. 

3)  DC  level-shifting.  Tests  the  effects  of  AC 
coupling  the  detector  outputs  on  image  reconstruction. 

The  routines  used  to  perform  these  operations  transformed  the 
detector  vector  disk  file  and  stored  the  conditioned  waveform  in 
a  local  memory  buffer.  The  codes  were  organized  as  signal 
conditioning  options  on  a  menu  at  the  sta**;  of  the  IMAGE  program. 

4.7.1  Waveform  shift 

The  detector  vectors  are  packed  into  a  linearly 
organized  file  Properly  unpacked,  the  detector  values  can  be 
treated  as  a  two-dimensional  array,  so  that  any  element  cl  can 
be  labeled  as  where  q  =  1  -  4  (quadrant  number)  and 

t  -  1  -  465  (time  sample  number). 

A  phase  shift  of  one  full  sample  can  be  accomplished 


by  the  cyclic  permutation: 
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d(q,c)  = 

a(q,t*l) 

t  1  -  464  positive 

phase 

d(q,4S5) 

=  ^;q. i ) 

s  h  i  f  t 

dtq.t)  = 

d  ( q ,  t- 1 ) 

t  =  2,465  negative 

phase 

d(q,l)  = 

d(q,465) 

shi  ft 

A  phase  shift  of  less  than  one  sample  is  clone  by 
sample  interpolation: 

d(q,t)  =  d  ( q ,  t +  $(d(q.t+l)  -  :l(q,t))  t  =  1,464 
d(q,465)  =  d(q,465)  +  s(d(q,I)  -  d(q,465)  ) 

positive  phase  shift 

d(q,t)  =  d(q,t)  +  s(d(q , t-1 }  -  d(q,t))  t  -  2,466 
d ( q » 1 )  =  d(q,l)  +  s(d(q,465)  -  d(q, 1) ) 

negative  phase  shift 

By  iterative  reconstruction  of  a  simulated  input, 
detector  vector  conditioned  by  various  phase  shifts,  the  sensitivity 
to  waveform  phase  was  determined.  A  shift  of  a  full  time  sample 
was  found  to  completely  destroy  the  reconstruct  ion .  The  largest 
phase  shift  that  can  be  used  without  unacceptable  image  degradation 
is  approximately  .1  time  sample.  Figure  4.13  shows  a  simulated 
horizontal  bar  phase  shifted  by  .1  time  sample.  The  total  power  in 
the  image  is  reduced  by  about  1,  e. 

4.7.2  Waveform  smoothing 

In  order  to  test  the  effect  of  waveform  smoothing  as 
might  be  caused  by  electronic  bandwidth  limitations,  a  simple 
digital  smoothing  filter  was  used: 


(>/ 

d  ( q ,  t )  =  q(d(q,t-l)  +  2d(q,t)  f  d(q,t  +  l)  )  t  +  2,464 
d  ( q » 1 )  =  14(d(q,465)  +  2d(q,l)  +  d(q,2)) 
d(q,465)  -  a(d(q,464)  +  2d(q  ,465)  +  d ( q  ,  1 ) ) 

The  result  of  smoothing  the  waveform  is  a  smoothing  of  the 
reconstructed  image.  Figure  4.14  shows  the  results  of  passing 
the  detector  vector  through  the  smoothing  filter  once,  while 
Figure  4.15  is  the  same  image  after  three  passes  of  the  waveform 
smoothing  filter.  The  major  effect  is  a  continued  loss  of  image 
resolution. 

4.7.3  Level  shifting 

The  AIT  reconstruction  algorithm  was  defined  with  the 
implicit  assumption  that  the  input  detector  samples  were  specified 
positive  definite  with  respect  to  an  absolute  zero  level.  From 
an  engineering  point  of  view,  this  implies  a  DC  coupled  detection 
system.  In  an  AC  coupled  system  the  detector  waveform  has  positive 
and  negative  value  extremes  such  that  the  RMS  value  of  the  waveform 
is  zero.  This  results  in  an  offset  or  negative  DC  pedestal  on  the 
waveform.  Figure  4.16  shows  the  results  of  level  shifting  the 
simulated  waveforms  for  the  horizontal  bar  negative  by  one-half 
the  peak  waveform  value.  The  reconstruction  of  the  bar  is  unaf¬ 
fected  but  artifacts  are  created  in  the  corners  of  the  reconstruction 
field. 

The  DC  offset  of  the  waveform  is  physically  equivalent 
to  superposing  the  target  bar  on  a  uniform  (negative)  intensity, 
including  intensity  outside  the  field  of  view  of  the  algorithm  as 
defined  by  the  pixel  locations  which  generated  the  forward  matrix. 

- - - - — — - — —4- —  mi'  . 


lift: 
I  J4 


This  power  outside  the  reconstruction  field  is  aliased  into  the 
field,  forming  the  corner  artifacts. 

In  order  to  employ  AC  coupling  in  the  AIT  imaging 
system,  a  simple  DC  restoration  step  was  used  to  re-reference 
the  waveform  to  zero  absolute  minimum.  This  allows  effective 
operation  of  the  reconstruction  in  either  AC  coupled  or  DC  coupled 
systems  exhibiting  offset  drift. 
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5 ■ 0  SOFTWARE  SIMULATION:  FRACKING 

The  gen  ralization  of  the  original  SCS  tracking  algorithm  to 
accommodate  the  more  complex  elliptical  nutation  pattern  necessitated 
the  development  of  new  software  to  perform  and  test  the  algorithm  and 
to  rest  the  tracking  concept. 

The  procedure  for  implementing  the  AIT  tracker  algorithm 
defined  in  Section  3  is  shown  in  the  flow  diagram  of  Figure  5.1. 

The  generation  of  al1  the  parameters  required  for  execution  of  the 
algorithm  was  coded,  so  that,  given  a  description  of  the  target 
object  and  the  nutation  minor  cycles,  the  mask  functions,  gains,  and 
offset  corrections  were  calculated. 

Initial  programs  written  to  test  the  tracking  algorithm 
created  detector  vectors  by  assuming  the  target  to  be  a  point 
source  with  a  given  displacement  relative  to  the  center  of  the 
field  of  view.  These  programs  would  also  determine  the  "allowable 
arcs"  to  be  used  fo>°  a  given  object  size  (even  though  the  target 
was  considered  a  point  source).  Using  the -^e  detector  vectors  and 
allowable  arc  definitions,  the  programs  would  also  calculate  the 
gain  and  offset  for  each  tracker  cycle.  Algorithm  consistency  was 
initially  verified  ty  using  the  gains  and  offsets  calculated  for  a 
given  object  displacement  in  determining  the  positions  of  point 
sources  with  different  displacements.  Thereby,  the  sensitivity 
of  the  algorithm  to  sets  of  displacements  measured  across  sets  of 
gains  and  offsets  was  examined. 
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The  15/31  collapsing  ellipse  notation  pattern  was  divided  into 
minor  cycles  as  shown  in  Figure  5.2  -  5.4. 
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5 • 1  Tracking  Performance 

Transfer  functions  for  an  AIT  tracker  algorithm  defined  using 
a  2  pixel  sized  generator  object  were  simulated.  Small  displacement 
response  was  investigated  using  various  test  object  sizes,  while 
large  displacement  response  was  investigated  over  the  entire  8x8 
field  using  single  pixel  test  objects. 

Figure  5.5  shows  the  tracker  transfer1  function  for  small 
displacements.  In  this  simulation,  a  roughly  two-pixel  diameter 
Gaussian  intensity  profile  spot  was  translated  along  the  x-axis 
in  the  range  -.75  pixels  to  +.75  pixels,  and  the  tracker  response 
was  calculated  for  several  different  minor  cycles.  The  results 
are  excellent.  The  residual  offset  spread  for  the  transfer  curves 
through  zero  is  less  than  .06  pixel  size,  or  less  than  .03  of  the 
2  pixel  spot,  diameter.  Furthermore,  the  average  of  the  ninor 
cycle  offsets  is  found  to  be  zero  to  within  the  accuracy  of  the 
simulation.  The  algorithm  displays  no  residual  bias  for  averaging 
times  of  one  full  pattern  length  or  longer. 

In  Figure  5.6  the  same  calculation  is  repeated  using  a  4  x  4 
pixel  square  test  object  with  Gaussian  smoothed  edges.  This  object 
is  larger  by  twice  than  the  lindt  allowed  by  the  definition  of  the 
algorithm.  Nevertheless,  operation  is  very  good,  with  about  .1 
pixel  offset  spread  and  *  10'  effective  gain  variation  in  the 
transfer  func. .on  for  different  minor  cycles.  Again  the  average 
offset  for1  the  minor  cycles  is  zero.  This  result  shows  that  tracker 
performance  degrades  gracefully  as  the  object  size  exceeds  teat  for 
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which  the  algorithm  coefficients  are  optimized,  which  is  a  very 
important  practical  feature. 

The  large  signal  response  curves  also  show  useful  performance 
features.  Figure  5.7  shows  the  tracker  response  to  a  one-pixel 
size  spot  translated  between  the  coordinates  x  =  +  3.5  pixels.  The 
path  is  displaced  to  offset  y  -  .  5  pixel.  Even  at  the  edges  of 
the  reconstruction  field  the  slope  change  are  not.  large. 

In  Figure  5.8,  the  object  is  again  translated  in  the  x-direction, 
but  with  a  displacement  of  2.5  pixels  along  the  y-axis.  Even  at 
the  extremes  of  the  displacement  the  curves  are  monotonic.  This 
property  holds  throughout  the  entire  field;  i.e.  there  are  no  slope 
reversals  that  would  result  in  any  system  instability.  The  algorithm 
under  all  conditions  responds  with  the  correct  sign  of  error  signal, 
cind  in  anv  closed-loop  application  would  drive  the  system  cowards 
null,  where  the  ultimate  response  is  excellent. 


WHMn 
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6 . 0  BREADBOARD  HARDWARE:  OPTICAL  SYSTEM 

As  part  of  the  AIT  program,  a  breadboard  optical  system  was 
designed  and  constructed  to  generate  nutated  detector  waveforms  from 
real  target  images.  The  images  were  produced  by  a  target  simulator 
in  which  target  size,  shape,  intensity,  position  and  rotation  could 
be  controlled.  Both  systems  were  designed  to  operate  at  either 
visible  or  infrared  wavelengths. 


6.1 


Optical  Layout 


Figure  6.1  is  a  diagram  of  the  AIT  optical  head  and  target 
simulator  system,  drawn  approximately  to  scale.  A  photograph  of  the 
set-up  is  shown  in  Fiqure  6.2. 

The  simulated  target  originates  as  a  target  mask,  shown  at  top 
center  in  Figure  6.1.  For  visible  operation,  the  target  could  be 
a  continuous-tone  35  mm  transparency  or  a  thin  chemically-etched 
copper  mask.  For  IR  operation,  only  the  copper  masks  would  be  used. 
The  target  is  back-illuminated  by  a  simple  tungsten  lamp  and  condenser 
arrangement;  in  the  (R,  the  tungsten  lamp  would  be  replaced  by  a 
black-body  source.  The  target  mask  itself  is  mounted  on  a  stepper- 
motor  driven  rotary  stage,  allowing  target  rotation  to  be  simulated. 

The  expanding  beam  from  the  target  projector  passes  through  a 
pair  of  galvanometer  mirror  deflectors,  shown  in  the  upper  left 
corner  of  the  figure.  The  deflectors  can  impose  electrically 
controlled  tip-tilt  on  the  beam,  inducing  apparent  motion  of  the 
target  object . 


AIT  TARGET  SIMULATOR  -  OPTICAL  BREADBOARD 


HONEYCOMB  BENCH 
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The  beam  is  folded  and  sent  to  an  off-axis  paraboloid,  which 
collimates  the  beam.  The  beam  is  then  reflected  from  a  flat  mirror 
mounted  or.  piezoelectric  actuators.  This  mirror  can  also  cause 
apparent  target  motion,  or  it  can  be  used  as  a  steering  mirror  in 
closed-loop  tracking  experiments . 

From  the  piezoelectric  tracker  mirror  the  beam  is  reflected 
to  another  off-axis  paraboloid,  which  bnngs  the  beam  to  a  focus  at 
the  input  of  the  optical  head  after  reflection  from  a  final  fold 
mi rror. 

The  optical  train  of  the  optical  head  must  provide  for  the 
following: 

1)  The  input  must  be  imaged  on  the  output  quad  cell 

detector. 

2)  The  nutation  pattern  must  be  impressed  on  the 

input  beam. 

3)  Means  for  adding  a  reference  i o  the  input  beam 
must  be  available. 

4)  Means  for  steering  the  input,  beam  from  the  inside 
of  the  instrument  in  an  accurate  grid  pattern,  with  an  eye  to 
generation  of  an  instrumental  back  matrix  tsec?  section  2.4.2)  is 
desi rable. 

The  refer  nee  sv  tern  was  desiunen  fur  : f{  use  of  the  hr  ea.lb.-a t  d , 
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constructed  to  inject  a  reference  beam  at  this  point. 

The  large  spherical  mirror,  shown  top  right  in  the  figure,  is 
the  primary  element  for  re-imaging  the  input  beam  on  the  detector. 

It  is  equipped  with  motor-driven  micrometer  actuators  and  a  precision 
tilt  sensor,  so  that  it  can  be  accurately  positioned  electronically 
for  generation  of  an  instrumental  back  matrix. 

In  this  system,  nutation  is  performed  by  a  set  of  galvanometer 
scan  mirrors,  which  separately  impose  the  x  and  y  nutation  deflections. 
In  order  that  the  pupil  be  imaged  on  each  mirror  without  intervening 
optical  elements,  an  anamorphic  field  lens  is  placed  at  the  input  of 
instrument.  This  causes  the  pupil  image  to  be  separated  to  form  to 
orthogonal  line  focii,  coinciding  with  the  rotation  axis  of  each 
nutator  mirror.  A  set  of  output  lenses  forms  the  nutated  target 
image  at  an  appropriate  magnification  on  the  silicon  quad  cel i 


detecto r. 


86 


6 . 2  Nutatj  on 

The  AIT  tracker-imager  requires  art  optical  scan  pattern  of 
sufficient  sampling  density  to  meet  the  spatial  Nyquist  criteria. 

The  rule-of-thumb  for  a  quad  cell -based  AIT  is  that  the  quad  cell 
axes  must  sample  each  pixel  at  four  distinct  positions  for  alias-free 
reconstruction  to  be  acnieved.  The  scan  pattern  used  in  Phase  I  was 
a  circular  spiral  which  contracts  linearly  in  time  from  its  maximum 
to  minimum  diameter,  where  the  latter  is  one-eighth  the  former. 

This  is  called  the  linear  spiral.  The  spiral  takes  eight  revolutions 
to  go  from  maximum  to  minimum  diameter  and  eight  more  to  return  to 
maximum,  for  „  +otal  of  sixteen  (16)  revolutions  per  image  cycle. 
During  each  revolution,  thirty-two  (32)  time  samples  are  taken  in 
each  quadrant,  for  a  total  of  32  x  16  x  4  =  2,048  data  points  per 
image. 

The  resulting  sampling  apttern  is  shown  in  figure  6.3a.  The 
spoke-like  appearance  of  the  pattern  results  from  the  use  of  uni¬ 
formly  spaced  time  samples.  For  the  same  reason,  the  sampling  density 
is  much  higher  near  the  center  of  the  pattern.  This  inefficient  use 
of  the  nutation  time  causes  the  exterior  pixels  to  he  marginally 
sampled  while  the  interior  is  over-sampled.  Figure  6. 3b  shows  the 
path  described  by  the  nutation  device  over  the  pixel  field. 

The  x-  and  y-drive  waveforms  for  this  pattern  are  shewn  as  the 
bottom  two  traces  of  Figure  6.3c.  The  signals  are  sine  and  cosine 
waves  multiplied  by  linear  shaped  envelope  functions,  which  are 
shown  for  x  and  y  as  the  first  and  second  traces  respectively. 
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For  the  experimental  work  done  in  Phase  I,  a  two-axis  piezoelectric 
deflector  driven  by  a  digitally  programmable  waveform  generator  was 
used  to  produce  this  nutation  pattern.  It  was  adequate  for  demon¬ 
stration  purposes  but  is  too  slow  (approximately  250  Hz  useful  fre¬ 
quency  range)  to  be  usable  in  a  field  instrument. 

3 

Previous  work  on  the  related  I  sensor  utilized  resonant 
galvanometer  scan  mirrors  to  provide  a  circular  nutation  pattern 
with  constant  amplitude  at  10  kHz  scan  rates.  Using  two  mirrors 
(one  per  axis)  in  phase  quadrature,  the  nutation  can  be  controlled 
to  better  than  1%  in  both  amplitude  and  phase. 

However,  the  resonant  scan  mirrors  are  designej  to  be  high-Q 
devices  (Q  is  greater  than  1000)  in  crder  to  achieve  the  required 
deflections  at  high  speeds.  At  10  kHz,  this  implies  a  antrol  band¬ 
width  of  less  than  100  Hz.  A  Fourier  analysis  of  the  triangle 
wave  envelope  function  of  Fiugre  6.3c  yields  frequency  components 
at  the  fundamental,  (10 7 16)  =  625  Hz,  with  higher  harmonic  terms 

at  frequencies  (625  t  625  x  n)  Hz.  n  odd,  whose  amplitude  decreases 
2 

as  1/n  .  This  means  that  the  mirrors  must  respond  tc  frequency 
components  several  kilohertz  above  and  below  the  10  kHz  center 
frequency  in  order  to  approximate  the  required  deflection  waveform!' 
(Figure  6.3c,  bottom).  This  is  precluded  by  the  high  mirror 
(}.  The  conclusion  is  that  a  pair  of  resonant  galvns  cannot 
generate  the  necessary  All  nutation  pattern. 
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6.2. 1  i our-mi rror  nutator 

A  solution  to  this  problem  is  shown  in  Figure  6.4.  Here 
the  scan  mirror  in  each  axis  is  replaced  by  a  closely  spaced  pair  of 
mirrors.  The  basic  ice a  is  to  drive  each  mirror  of  a  one-axis  pair 
with  a  different,  constant,  frequency.  The  amplitude  modulated 
deflection  waveform  for  that  axis  is  then  given  by  the  beat  frequency 
modulation  between  tho  two  mirror  drives.  One  way  to  view  this 
method  is  to  consider  a  drive  waveform  containing  only  two  Fourier 
components  (instead  of  many,  as  with  triangle  wave  modulation). 
Instead  of  attempting  to  force  a  single  mirror  to  respond  to  this 
waveform,  the  components  are  separately  applied  to  two  mirrors  and  th 
superposition  obtained  optically. 

The  advantage  of  this  system  is  that  all  mirrors  can  be 
run  at  a  constant  phase  and  amplitude.  The  techniques  for  accom¬ 
plishing  this  are  already  demonstrated. 

6.2.2  Four-mirror  nutation  patterns 

Two  classes  of  nutation  patterns  can  be  obtained  using 
the  four-mirror  system,  referred  to  as  the  cosine  spiral  and  the 
collapsing  ellipse.  They  can  be  defined  by  the  drive  waveforms  for 


each  mirror: 
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Pour-mirror  nutator  geometry. 


FIGURE  6.4 
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Cosine  spiral  nutation: 

Xj(t) 

1 

’  W2  ' 

^sin2ir  j 

>  +  <V2S>)  lS 

X2(t) 

*  W2  j 

fsin2tt[ 

L.  " 

>N  -  <V2S>J  ’1 

y,(t) 

•  W2 1 

cos2“ft  | 

;N  +  (v2si3  ^ 

y2(t) 

■  “max/2  \ 

!  cos2T[V 

ti  -  (v«U 1  \ 

x-|(t),  x2(t)  ~  waveform  drives  for  x-axis  mirrors 

y-j(t),  y2( t)  =  waveform  drives  for  y-axis  mirrors 

Dm.„  =  maximum  diameter 
max 

fN  =  nutation  frequency 

S  =  number  of  spirals  per  pattern  (S/2  in,  S/2  out) 


The  composite  x  and  y  scan  amplitudes  will  be  (  by  a 
simple  trigonometric  substitution:): 


“spiral  '  °™x!in^fNt 


y  .  ,  =  D  cos2«fwt  cos2irfw  /ct 

Jspiral  rax  N  .  N/S 
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Collapsing  eMipse  nutation: 


X-,  (»t ) 
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fN 
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Although  the  cosine  spiral  and  ellipse  waveforms  differ 
only  in  the  phase  of  the  y-axis  envelopes,  the  resulting  patterns 
are  qualitatively  very  different. 

Figure  6.5a  shows  the  sampling  density  for  a  cosine  spiral 
with  S  =  15  revolutions  per  pattern  and  31  time  samples  per  revo¬ 
lution.  The  nutation  path  is  that  shown  in  Figure  6.5b  while  6.5c 
gives  the  x-and  y-  envelope  functions  and  drive  waveforms.  Although 
the  edge  sample  density  is  better  than  for  the  linear  spiral  of 
Figure  6.3a,  the  sampling  at  intermediate  radius  is  somewnat 
sparse,  while  the  central  pixels  are  still  over-sampled.  The  cosine 
spiral  would  still  have  performance  comparable  to  that  of  the  linear 
spiral,  and  is  realizable.  It  remains  a  poor  match  for  a  square 
reconstruction  field,  and  so  the  following  pattern,  the  collapsing 
ellipse,  was  used  instead. 
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Figure  6.6a  shows  the  sampling  pattern  for  the  collapsing 
ellipse,  for  S  =  15  revolutions  per  pattern  and  31  samples  per 
revolution.  The  square  fill  shape,  typical  of  Lissajous  figures,  is 
an  ideal  match  to  a  square  pixel  field,  and  the  sampling  density  is 
now  weighted  toward  the  field  edges.  Figure  6.7  matches  this 
sampling  pattern  to  an  8  x  8  pixel  field.  It  is  seen  that  the  rule- 
of-thumb  of  four  distinct  samples  per  pixel  is  obtained  for  every 
pixel  in  the  field.  Thus,  the  collapsing  ellipse  with  these  scan 
parameters  yields  a  nearly  ideal  nutation  pattern  for  image  recon¬ 
struction.  Inspection  of  the  individual  mirror  drive  waveforms 
shows  also  that  the  pattern  is  generated  using  only  sine  components, 
eliminating  the  need  to  maintain  a  quadrature  phase  reference,  thus 
simplifying  the  nutation  controller.  A  set  of  galvanometers  was 
designed  and  constructed  vhich  could  operate  at  a  basic  nutation 
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6 . 3  Interface  Electronics 


In  order  to  test  the  imaging  capabilities  of  the  AIT 
breadboard,  a  set  of  special  interface  electronics  was  constructed 
to  allow  the  breadboard  to  operate  directly  with  the  DG  Nova  computer. 
This  enabled  development  and  debugging  of  the  optical  subsystems  in 
parallel  with  work  on  the  PMP. 

The  interface  system  is  diagrammed  in  Figure  6.8.  The  four 
silicon  quad  cell  outputs  are  buffered  by  a  set  of  preamplifiers, 
and  the  analog  waveforms  sent  to  a  set  of  analog  switched  integrators, 
The  integrators  performed  the  front  end  signal  integration  and  sampling 
required  by  the  AIT  algorithm,  under  the  control  of  timing  signals 
from  the  DG  Nova.  The  output  sampled  waveforms  were  digitized  by  the 
Nova  DG/DAC  and  recorded  on  disk. 

In  addition  to  data  conditioning,  the  interface  package 
controlled  the  AIT  slow-speed  galvo  set.  For  debugging  and 
experimentation  without  the  PMP,  a  linear  response  pair  of  galvos, 
with  sufficiently  wideband,  albeit  slow,  response  was  used,  since 
the  Nova  DG/DAC  was  not  fast  enough  to  accommodate  the  7980  Hz  rate 
of  the  fast  A  mirror  galvo  set.  The  slow  galvos  were  driven  by  the 
compound  waveforms  of  Figure  &.6c,  bottom.  1  he  waveforms  were 
generated  on  an  Apple  computer  and  downloaded  to  the  interface  wave¬ 
form  generator  through  an  RS-232  interface.  The  waveforms  were  ap¬ 
plied  to  a  set  of  driver  amplifiers  which  controlled  the  linear 
gal vos . 
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7 •  0  BREADBOARD  >1ARPWARE :  ELECTRONIC  PROCESSOR 

The  complex  AIT  imaging  and  tracking  algorithms  require  advanced, 
highly  specialized  digital  hardware  to  realize  real-time  operation. 

To  this  end,  a  custom  digital  signal  processor,  the  Programmable 
Microcoded  Processor  (PMP),  was  designed  and  built, 

7 • 1  Processor  Requirements 

The  processor  unit  for  the  AIT  system  must  be  capable  of  per¬ 
forming  imaging,  tracking,  nutation  control,  system  control,  data 
output  and  other  tasks  in  real  time.  This  involves  accessing  of 
detector  outputs,  and  sampling,  digitizing  and  buffering  and 
organizing  the  results  into  two  data  sets:  the  major  cycle  data 
vectors  which  are  1860  samples  long  and  must  be  processed  for  image 
reconstruction,  and  the  minor  cycle  sub-vectors  124  samples  long, 
which  must  be  converted  to  tracker  outputs.  Tracking  must  ultimately 
occur  at  3  kHz  update  rate,  while  the  imaging  goal  is  30  frames/sec. 

7 . 2  PMP_  An c h  i  tec  tu re 

The  very  high  data  throughput,  rate  required  of  the  system  pro¬ 
cessor  precludes  use  of  any  form  of  standard  minicomputer  architec¬ 
ture.  Tnese  general  purpose  machines  require  decoding  of  each 
instruction  into  the  hardware  commands;  the  decoding  step  takes 
several  machine  cycles,  greatly  diminishing  sneed.  I ven  advani ed 
array  processors  are  encumbered  Ipy  use  of  high- accuracy ,  float  ing 
point,  arithmetic  unneeded  here,  as  well  as  cum!  'some  I/O  structures, 
again  reducing  thruuqhpu  1  he  alternative  is  a  dedicated  at  chi  t.ec  furr 
with  no  instruction  decoding,  fixed-point  arithmetic  and  optimized 
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I/O,  which  meets  performance  requi remen ts ,  and  is  also  at  least  a 
factor  of  ten  less  expensive  than  the  closest  commercially  available 
processor  configuration. 

The  PMP  is  a  purely  microcode  driven  digital  signal  processor 
with  16  bit  data  and  address  busses  and  16  bit  basic  fixed  point 
calculation  accuracy.  Based  on  the  AMD  2903  family  of  bit-slice 
components,  the  system  has  a  100  nS  cycle  time  and  employs  a  96 
bit  fully  horizontal  microcode  'word  size. 

A  simplified  diagram  of  the  PMP  architecture  is  shown  in  Figure 
7.1.  It  is  apparent  that  the  system  is  complex;  the  speed  and 
flexibility  of  the  processor  have  been  obtained  at  the  expense  of  a 
complicated  architecture  which  is  challenging  to  program.  The  effort 
has  nevertheless  been  highly  successful,  resulting  in  a  very  reliable 
system  without  a  single  analog  adjustment. 

The  PMP  can  be  approached  in  terms  of  functional  sub-units  or 
processing  resources.  The  extensive  use  of  tri -state  bus  switches 
allows  different  groups  of  hardware  components  to  function  quasi¬ 
independent  ly  on  different  tasks,  with  the  boundaries  between  groups 
and  their  communication  channels  variable  dynamically  under  micro¬ 
program  control. 

At  center  left  in  the  diagram  are  the  Am  2903  and  2902  units 
which  form  the  system  ALU.  They  are  loaded  and  operated  by  the 
component  group  at  top  center,  the  Am?904  and  Am?910,  the  96  bit 
wide  microcode  memory  and  trip  microcode  pipeline  register.  fhese 
devices  form  a  powerful  control  unit  for  directing  data  flow, 
memory  access,  external  device  operation  and  internal  computation 
with  a  significant  degree  of  parallelism.  Ihe  efficient  control  of 
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computational  flow  is  enhanced  by  the  group  of  four  memory  units 
shown  below  the  Am2910  group:  the  Bus  Address  Register  (BAR)  and  its 
Direct  Bus  Access  (DBA)  unit,  and  the  Memory  Address  Register  (MAR) 
and  Direct  Memory  Access  unit.  These  memories  contain  the  maps  and 
tables  needed  to  fetch  and  route  data  in  the  order  required  from 
the  main  data  memory  without  use  of  any  time-consuming  address 
cal cul ations . 

The  main  data  memories  are  located  at  lower  left.  They  consist 
of  a  pair  of  4K  by  16  bit  blocks  with  separate  access  ports  and  address 
controls.  This  allows  one  bank  to  inout  data  from  some  external 
source  (such  as  the  digitized  output  af  a  detector  set)  while  the 
other  is  read  out  for  processing,  thus  implementing  a  "ping-pong1" 
buffer. 

In  addition  tc  these  major  PMP  hardware  units,  the  AIT  version 
of  the  PMP  contains  several  other  sub- sys terns ,  shown  at  right  in 
the  figure.  The  coefficient  memory  contains  128K  by  16  bit  words 
of  storage  to  accommodate  the  back  matrix  for  image  reconstruction. 
Inputs  from  the  detector  preamps  are  integrated  and  digitized  in 
the  discrete  samples  needed  by  the  imaging  and  tracking  algorithms 
in  a  detector  A/D  front  end  with  digitally  controlled  AGC,  shown  at. 
lower  right. 

Generation  of  up  to  four  separate  nutation  drives  and  timing 
of  all  nutation-synchronized  operations  is  done  by  the  Timing  and 
Image  memories,  located  at  right-center.  Analog  outputs  of  essentially 
any  internal  values  or  results,  such  as  x  and  y  tracker  outputs,  are 
facilitated  by  a  bank  of  eight  R-bit  D/A  converters. 
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finally,  communication  v/i  t.h  a  host  or  supervisory  computer  must 
be  considered.  Since  the  PMP  is  a  dedicated  digital  system  deliberately 
containing  no  high-level  instruction  capability,  operator  access, 
system  downloading,  and  other  such  tasks  must  be  accomplished  through 
some  external  general  purpose  computer.  This  communication  proceeds 
through  a  standard  RS-232  port,  on  the  data  bus.  Through  this  port 
any  and  all  memory  locations  and  registers  (which  are  equivalent  in 
the  PMP  addressing  scheme)  can  be  loaded  or  accessed;  in  particular, 
the  IK  by  96  bit  microcode  memory  is  downloaded  through  this  port. 

Major  support  of  the  PMP  was  done  through  the  Data  General  Nova, 
on  which  all  microcode  was  developed  using  a  specially  configured 
rneta-assp' ■'  r.  For  system  diagnostic  purposes,  Pascal  code  for  an 
Apple  computer  was  also  written  whicn  allowed  it  to  act  as  an 
intelligent  terminal  for  the  PMP.  These  support  systems  helped  "”ake 
the  challenging  task  of  programming  the  PMP  more  manageable  and 


eff  cient. 
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A IT  Mi  crop rogram  Stru e  tu  re 
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The  AIT  microcode  procedures  are  divided  into  a  set  of  modules 
contained  within  a  simple  execution  loop.  The  modules  are  used 
to  sequentially  initialize  each  of  the  PMP  memories  with  the  algorithm 
parameters  (tracker  mask  function  and  offsets,  nutation  waveforms, 
etc.)  required  for  a  particular  AIT  configuration.  As  a  first  step, 
a  16-word  bootstrap  loader  is  entered  into  the  top  of  microcode 
memory,  which  facilitates  loading  of  the  remaining  micro-program. 

The  operation  of  the  AIT  microprogram  is  outlined  in  Table  7.1, 
while  the  flow  of  program  control  is  diagrammed  in  Figure  7.2.  The 
action  of  the  microprogram  proceeds  as  follows. 

After  the  PMP  has  been  downloaded  with  the  All  microcode,  the 
microprogram  will  enter  a  loop  waiting  for  a  word  to  be  transmitted 
from  the  host  over  the  RS-232  line,.  When  the  word  is  received,  it 
will  be  interpreted  as  a  command  by  the  microporgram,  which  will 
branch  to  a  specific  module  using  the  command  code  as  an  offset  into 
a  dispatch  table.  After  the  command  has  been  executed,  control  will 
return  to  the  wait  loop  (except  for  the  "RUN"  command)  where  the 
mircoprogram  will  wait  for  the  next  command.  Thus,  the  AIT  micro¬ 
program  contains  no  executive  (only  a  central  dispatcher  and  several 
modules),  and  is  totally  driven  by  the  host.  This  allows  the  PMP  to 
operate  in  a  very  general  environment  with  the  PMP  run-time  configura¬ 
tion  dictated  by  the  hosi  to  allow  running  different  variations  of 
the  same  experiment. 
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AIT  SCENARIO 

.  DOWNLOAD  TiiAGE  MEMORY.  Contains  digital  codes  to 
be  convened  to  analog  signals  which  drive  the  two 
nutation  mirrors. 

.  DOWNLOAD  TIMING  MEMORY.  Contains  digital  timing 
and  synchronizing  information  used  to  integrate  and 
digitize  QUAD  CELL  outputs,  ana  synchronize  micro¬ 
program  to  the  data  acquisition  process. 

•  DOWNLOAD  DBA  MAP.  Contains  codes  used  to  address  and 
read  the  integrating  analog  to  digital  converters. 

•  DOWNLOAD  DAR  MAP.  Contains  address  pointers  used 
to  fetch  the  digital  samples  and  pre-stored  coeffi¬ 
cients  to  allow  high  speed  algorithm  execution. 

•  DOWNLOAD  MAIN  MEMORY.  Contains  algorithm  constants. 

•  DOWNLOAD  COEFFICIENT  MATRIX  MEMORY.  Contains 
Imaging  Transformation  Matrix  (Back  Matrix). 

•  RUN. 

1.  Initialize  and  control  data  acquisition  hardware. 

2.  Execute  Tracking  algorithm. 

3.  Execute  Imaging  algorithm. 


TABLE  7.1 


The  PMP  contains  three  autonomous  dedicated  hardware  resources 
that  need  to  be  initialized.  One  is  the  Image  Memory  which  contains 
the  digitized  nutation  control  waveforms  used  to  drive  the  mirrors. 

The  host  computer  will  generate  these  waveforms  and  send  the  sequence 
of  numbers  to  the  PMP,  which  will,  in  turn,  load  the  sequence 
into  the  memory.  In  addition,  the  microprogram  will  initialize  the 
Image  Memory  control  hardware  that  is  used  to  determine  the  nutation 
frequency.  Once  initialized,  the  Image  Memory  machine  is  free-rurining, 
and  the  waveform  sequence  will  be  repetitively  fetched  and  converted 
to  analog. 

The  second  resource  is  the  Timing  Memory  which  contains  timing 
information  used  to  synchronize  all  the  elements  in  the  system. 

One  bit  is  used  to  command  the  A/Ds  to  convert  and  another  bit  is 
used  to  synchronize  the  microprogram  to  specific  points  in  the 
nutation  cycle  (i.e.  tracker  cycle  endpoints).  Again  the  Timing 
Memory  is  loaded  by  the  microprogram  with  a  timing  map  generated  by 
the  host;  and  its  control  logic  initiated  so  that  these  two  memories 
are  running  in  lock  step  with  a  programmable  phase  delay  between  them. 
Note  that  the  Timing  Memory  allows  the  microprogram  to  collect  and 
operate  on  the  sampled  data  in  real-time  and,  thus,  optimize  the 
tracking  response  time. 

The  third  resource  is  the  Direct  Bus  Access  (DBA)  machine  used 
to  fetch  the  four  integrated-digi tized  QUAD  CELL  outputs,  under  control 
of  a  Timing  Memory  bit.  The  detector's  outputs  (at  specific  points 
in  the  nutation  cycle,  with  minimum  latency)  are  stored  into  a 
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buffer,  in  a  fixed  order,  to  be  processed  later  by  the  microprogram. 

As  in  the  previous  two  memories,  this  dedicated  hardware  resource 
must  be  programmed  by  the  microcode. 

In  order  to  facilitate  high-speed  algorithm  execution,  another 
hardware  resource  exists  called  the  Direct  Access  Register  (DAR)  map. 
Essentially,  this  machine  maps  a  count  sequence  into  a  random 
sequence.  In  other  words,  a  sequence  generated  by  a  counter  accesses 
a  RAM  (containing  an  address  pointer  list)  which  in  turn  points  to  a 
fixed  storage  location  in  the  main  buffer  where  the  data  samples  are 
stored;  and,  thus,  obviating  any  effective  address  calculations. 

Also  stored  in  the  main  data  base  are  various  algorithm  parameters 
such  as  the  number  of  segments  in  a  tracker  cycle,  gain  coefficients, 
and  displacement,  offsets.  Whenever  the  microprogram  needs  to  access 
a  data  sample  or  a  parameter,  the  microinstruction  will  issue  the 
micro-orders  to  access  the  buffer  using  the  DAR  map  and  advance  the 
DAR  counter. 

The  only  unique  command  is  the  "RUN"  command  where  the  micro¬ 
program  executes  the  tracking  and  imaging  algorithms.  Note  that  these 
two  algorithms  are  independent  in  that  they  operate  on  different 
data  sets  and  under  different  constraints,  but  with  some  resource 
sharing.  Thus,  the  environment  is  similar  to  mul tiprogramming  with 
the  tracking  routine  having  higher  priority. 

Essentially  the  microprogram  will  periodically  test  a  bit  in 
Timing  Memory  signifying  that  all  the  segments  of  a  Tracker  Cycle 
have  been  integrated,  digitized,  and  stored  in  main  memory,  and  ready 
to  be  processed.  When  this  condition  occurs,  the  microprogram  will 


109 

branch  to  the  tracking  module.  The  microcode  first  fetches  the 
parameter  that  contains  the  number  of  segments  in  the  current  Tracker 
Cycle  and  loads  this  number  into  the  microsequencer's  internal  loop 
counter.  Next,  using  the  DAR  Map  and  the  loop  counter,  the  segments 
are  accumulated  for  each  of  the  detector's  four  quadrants.  Then 
the  routine  computes  the  "numerator"  and  "denominator",  performs 
the  division,  multiplies  the  ratio  by  the  gain  coefficient,  and  then 
adds  in  the  offset.  Finally  the  displacement  is  checked  for  over¬ 
flow  and  clamped  to  full  scale  if  necessary. 

Whenever  the  microprogram  is  not  executing  the  tracking  algo¬ 
rithm,  and  a  full  frame  of  "fresh"  data  has  been  collected,  the 
microcode  will  execute  the  imaging  algorithm.  Note  that  main  memory 
is  really  divided  into  two  separate  ping-pong  buffers  so  that  where 
it  is  time  to  process  a  new  image  frame,  the  buffers  are  swung  for 
one  major  cycle  in  order  to  freeze  one  frame's  worth  of  data.  On 
the  next  cycle  the  buffers  are  swung  back  so  now  the  imaging  data  can 
be  accessed  without  contention  over  several  major  cycles. 

The  imaging  algorithm  basically  consists  of  computing  the  inner 
product  of  the  detector  vector  (1860  sample  points),  times  the  Imaging 
Transformation  Coefficient  Matrix  (64  x  I860)  to  yield  the  pixel 
vector  (64  points).  First  the  detector  vector  element  is  fetched 
from  memory  (using  the  DMA  address  counter)  and  the  coefficient  is 
fetched  from  the  128K  C-matrix  memory,  then  multiplied  together 
(16  x  16  signed  multiply):,  and  finally  added  to  the  previous  product 
(32-bit  accumulation).  After  the  entire  column  has  been  multiplied, 
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the  pixel  value  is  sent  to  the  host  over  the  RS-232  line.  Note  that 
the  DC  bias  must  first  be  removed  from  the  detector  vector  before 
multiplication.  The  DC  bias  value  for  each  quadrant  is  specified  by 
choosing  those  four  elements  in  the  detector  vector  for  which  it  is 
guaranteed  no  power  fell  on  the  quadrant. 
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7.4  Processor  __Performaiice 

The  PMP  was  rtiicrocoded  to  perform  both  the  imaging  and  tracking 
funcitons  and  was  tested  using  simulated  data  sets  downloaded  from 
the  DG  Nova.  The  processor  performed  botn  tasks  flawlessly. 

As  a  part  of  system  operation,  the  PMP  had  to  generate  the 
15/31  nutation  pattern  drives  for  the  nutation  mirrors.  Figure  7.3 
shows  an  x-y  oscillogram  of  the  waveform  outs  of  the  drives  for  a 
two-mirror  system.  Amplitude  modulation  of  the  oscilloscope  was 
controlled  by  the  sample  timing  pulse.  The  result  is  an  accurate 
rendition  of  the  15/31  collapsing  ellipse  pattern.  This  verifies 
autonomous  operation  of  the  Image  and  Timing  memories  and  associated 
nutation-synchroni zed  ci rcu i try. 

To  test  the  full  range  of  tracker  operation,  an  entire  forward 
matrix  was  loaded  into  coefficient  memory,  where  the  simulated 
pixel  vectors  for  each  of  the  64  pixel  locations  could  be-  accessed  by 
the  PMP  in  the  same  way  as  real-time  data.  The  PMP  was  programmed 
to  step  through  each  pixel  in  raster  fashion,  convert  the  data  to 
spot  location  using  the  AIT  tracker  algorithm,  and  output  the  x,y 
results  througli  a  pair  of  the  general  purpose  D/A  converters. 

The  separate  X  and  Y  analog  outputs  are  shewn  in  figure  7.4. 

In  7.4a  the  full  scan  is  shown,  with  the  upper  trace  indicating  the 
x-output  and  the  lower  the  y-output.  The  expected  staircase  wave¬ 
forms  are  clearly  in  evidence,  reproducing  the  pixel  locations  to 
an  accuracy  consistent  with  the  tracker  simulation  results.  Figure 
7.4b  shows  one  row  of  eight  pixels  (constant  y  value),  revealing  the 
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departures  from  ideal  output  at  the  field  extremes  predicted  by 
the  simulations. 

Figure  7.5  shows  the  results  for  the  entire  field,  with  the 
x  output  plotted  against  y.  The  grid  pattern  associated  with  ideal 
operation  is  found  in  the  center  of  the  field,  with  loss  of  accuracy 
occurring  in  the  field  corners.  Again  this  matches  the  predictions 
of  the  simulations  and  shows  the  excellent  performance  of  both  the 
processor  and  the  algorithm. 

To  test  the  imaging  microcode,  the  back  matrix  was  loaded  into 
coefficient  memory,  and  simulated  detector  vectors  were  loaded  into 
data  memory  through  the  RS-232  port.  Two  vectors  could  be  loaded 
simultaneously,  one  in  each  buffer,  and  alternately  processed  using 
the  image  reconstruction  algorithm. 

The  results  were  found  to  be  excellent,  with  reconstruction 
fidelity  using  the  PMP’s  16  bit  fixed  point  arithmetic  essentially 
indistinguishable  from  the  DG  Nova  floating-point  results.  Image 
reconstruction  using  the  PMP  alone  took  approximately  2  seconds;  down 
from  the  16  required  by  the  DG  Nova  with  optimized  code.  Again, 
perfect  operation  of  the  PMP  was  demonstrated.  An  increase  of  frame 
speed  beyond  this  point  requires  use  of  the  special  hardware  multi¬ 
plier  unit,  the  Vector  Matrix  Multiplier  (VMM). 


X-Y  PLOT  OF  I  WO 
NUfATI ON  DRIVE  SIGNALS 
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7.5  Vector  Ma t r i x  Mu  1  tipi i or 

The  Vector  Matrix  Multiplier,  or  VMM  system,  is  made  up  of 
two  units.  One  unit,  called  the  Vector  Processor,  consists  of  a  dual 
vector  buffer,  a  vector  maf-mx  multiplier,  and  the  VMM  system  con¬ 
troller.  The  other  unit,  sailed  the  Coefficient  Memory,  consists 
of  an  expandable  coefficient  matrix  memory  system,  as  discussed  in 
Section  7.2.  Together,  these  units  perform  a  vector  matrix  multi¬ 
plication  for  a  single  row  vector  with  a  coefficient  matrix  with 
the  proper  dimensions.  For  the  AIT  study,  the  row  vector  has 
1  x  I860  elements  and  the  coefficient  matrix  has  1860  x  64  elements. 
After  a  vector  matrix  multiplication  is  performed,  a  new  vector  of 
64  elements  is  produced. 

VMM  Operating  Description 

Figure  7.6  describes  the  VMM  architecture.  It  is  divided  into 
a  Vector  Processor  and  Coefficient  Memory;  of  which  a  maximum  of  eight 
memory  sections  or  512K  words  of  memory  can  be  used.  Both  the 
Vector  Processor  and  Coefficient  Memory  have  been  designed  in  a  pipe¬ 
line  architecture  in  order  to  speed  the  vector  processing  through¬ 
put  rate.  The  maximum  through-put  rate  was  designed  for  10  MHz,  but 
at  the  present,  time  the  through-put  rate  is  l  imited  to  6  MHz  due  to 
the  multiplier  accumulator  device  being  used  (TRW  MAC  I010J).  As 
faster  multiplier  accumulator  devices  are  available,  the  maximum 
through  pul  rate  of  10  Mil/  can  be  realized. 
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In  operation,  the  VMM  vector  buffers  will  accept  data  from  the  PMP 
external  bus  (x  bus)  and  perform  the  multiplication  and  accumulation 
using  the  TRW  multiplier.  The  VMM  operates  unde;  its  own  clock, 
so  the  ping-pong  buffering,  with  the  two  sections  controlled  by  sepa¬ 
rate  address  generators,  acts  as  a  synchronous  interface  between  the 
PMP  and  VMM.  The  coefficients  are  fetched  from  C-momory  via  a 
separate  internal  bus  (in  contrast  to  the  way  C-memory  is  used  for 
PMP  image  processing)  to  prevent  tying  up  the  PMP  Data  bus  by 
coefficient  transfers.  The  resultant  image  vectors  are  sent  back 
to  the  PMP  over  the  PMP  data  bus. 

The  target  performance  of  the  VMM  is  shown  in  Table  7.2.  The 
hardware  is  designed  to  exceed  the  30  Hz  frame  rate  goal,  allowing 
a  healthy  operating  margin.  Although  completely  designed  arid 
fabricated,  the  VMM  was  not  debugged  and  integrated  with  the  PMP  by 
the  end  of  the  program,  and  so  full  30  frame/sec  processing  rate 
could  not  be  demonstrated. 


™^"^T'^^Wi^WWIiiW^PW,^^^  -^Tww^.n r.n' ^ — T~.r-^.-r  ..„TOr,!w-.,.- 


119 


•  Time  required  to  process  a  vector  with  one  column  of  the 
coefficient  matrix: 

1860  VMM  clock  cycles  (number  elements  in  vector) 

8  VMM  clock  cycles  (pipeline  feed  through  delay) 

TOTAL  1868  VMM  clock  cycle  (VMMCLK  -  6  MHz) 

TOTAL  TIME!  =  312  microseconds 

»  Time  required  to  read  one  multiplication  process  out  of  the 
VMM  into  the  PMP: 

10  PMP  clock  cycle 

The  PMP  clock  period  must  be  1/2  of  the  VMM  clock  period  to 
read  the  processed  result  (3  MHz) 

Therefore,  the  time  to  read  the  results  and  reinitialize 
another  multiplication  process  is: 

.33  microseconds  x  10  clock  cycles 

3.3  microseconds/readout 

•  Total  time  to  process  the  entire  matrix. 

(312  +  3.3)  64  =  20.2  milliseconds 

•  Conclusion: 

The  VMM  can  process  49.5  vector  matrix  multiplications 
per  second. 


TABLE:  7.2  VMM  SYSTEM  TIMING  REQUIREMENT  FOR  AIT 
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8 • 0  EXPERIMENTAL  SOFTWARE 

A  large  body  of  software  was  written  to  enable  the  Nova  to 
support  all  of  the  experimental  data  collection,  breadboard  develop¬ 
ment,  and  processor  development  tasks  in  the  program,  apart  from 
the  software  designed  for  algorithm  simulation.  In  the  sections 
below,  the  capabilities  of  the  codes  written  to  support  the  optical 
benchwork,  PMP  development,  the  Genisco  color  display,  and  general 
system  work  are  outlined.  The  descriptions  are  by  no  means  complete, 
but  rather  indicate  the  general  intent  and  emphasis  of  the  software 
effort.  A  listing  arid  brief  synopsis  of  all  the  programs,  sub¬ 
routines  and  utilities  written  for  the  AIT  program  is  given  in 
Appendix  I. 
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8 • 1  PMP  Support 

Programs  were  written  to  provide  four  types  of  support  to 
the  effort  to  integrate  the  AOA-designed  programmable  microcoded 
processor  (PMP)  into  the  AIT  project: 

1)  Programs  to  download  the  "bootstrap"  microcode  and 
the  microcoded  program  into  the  control  store  of  the  PMP. 

2)  Programs  to  interact  with  the  microcoded  diagnostic 
routines  to  debug  the  PMP  hardware.  There  is  a  "heirarchy"  of 
these  program  sets.  The  testing  proceded  from  verifying  the  inte¬ 
grity  of  the  most  vital  and  basic  hardware  resources  to  trying  to 
induce  subtle  and  data-dependent  hardware  failures. 

3)  Programs  to  simulate  the  microcoded  processes  using 
algorithms  that  mimic  these  processes.  Intermediate  results  could 
then  be  generated  and  displayed  and  used  as  an  aid  in  debugging 
the  microcode. 

4)  Programs  to  invoke  particular  microcode  modules  and 
to  process  and  display  the  results  calculated  from  the  data 
received  from  the  PMP. 

All  communication  between  the  PMP  and  the  Nova  is  carried  over 
a  standard  three-wire  RS-232  line.  A  very  simple  standard  protocol 
was  developed  for  communication  between  the  PMP  and  the  Nova.  This 
protocol  was  sufficient  for  all  the  modes  of  PMP  support.  Under 
this  protocol,  the  Nova  downloads  the  bootstrap  and  microcode 
program  to  the  PMP,  indicating  which  microcode  procedure  is  to  be 
carried  out.  If  more  data  is  required  from  the  Nova  for  the  PMP 
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to  carry  out  this  task,  more  information  is  sent  in  this  form:  a 
number  indicating  the  number  of  words  in  the  further  communication 
is  sent  to  the  PMP,  followed  by  the  data  words  themselves.  Finally, 
a  checksum  of  all  the  previous  words  is  sent  to  the  PMP.  The  PMP, 
in  order  to  communicate  with  the  slower  Nova,  awaits  a  one  byte 
request  from  the  Nova  before  sending  a  stream  of  data  (of  a  well 
defined  type  and  extent)  to  the  Nova. 

This  protocol  allowed  the  microcode  programs  to  be  built 
gradually  in  modules,  which  could  be  individually  called  and  tested 
by  the  Nova.  The  modular  structure  of  the  PMP  code  permitted 
greater  efficiency  in  coding,  and  also  permitted  greater  flexibility 
in  revising  arid  using  the  code.  Functional  microcode  modules  could 
be  selected  from  a  library  of  modules  and  matched  to  a  particular 
task  with  minimal  overhead. 

At  the  beginning  of  the  effort  to  achieve  PMP-Nova  integration, 
a  communications  problem  was  discovered.  The  Nova  RDOS  operating 
system  could  only  recognize  numbers  corresponding  to  ASCII  "control¬ 
s''  and  "control -Q"  as  console  control  characters,  even  when  the 
user  program  was  employing  a  Fortran  "read  binary"  command.  An 
assembly  language  routine  using  the  RDOS  system  command  ".RDS" 

(read  sequential)  was  written,  but  these  ASCII  characters  were  still 
not  recognized  as  data.  Both  input  techniques  resulted  in  these 
characters  being  intercepted  by  the  operating  system  and  thus, 
being  lost  to  the  calling  program. 
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Faulty  RDOS  documentation  (revealed  by  a  later  Data  General 
documentation  update)  prevented  a  simple  solution  to  this  problem 
From  becoming  immediately  apparent.  A  stopgap  measure  was  developed 
and  exploited  for  the  remainder  of  the  AIT  project.  This  solution 
consisted  of  writing  an  assembly  language  routine  to  directly 
interrogate  the  Universal  Line  Multiplexor  (ULM)  board,  the  hardware 
intermediary  for  coiTsnurti cations  with  the  PMP  and  the  Dasher  printer. 
This  ULM  interrogation  subroutine  was  used  in  all  the  Nova  programs 
communicating  with  the  PMP.  With  a  simple  revision  of  these  pro¬ 
grams,  ordinary  Fortran  "read  binaries"  could  now  be  used.  Many  of 
the  utilities  listed  in  Appendix  I  were  written  to  diagnose  and 
correct  this  problem. 
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8 *  2  Bench  Support 

Support  for  experimental  data  collection  was  provided  by 
programs  to  control  the  operation  of  the  electronic  devices  used  to 
drive  the  nutation  mirrors  and  to  sample  the  detector  outputs  in 
synchrony  with  the  sample  intervals  required  by  the  tracking  and 
imaging  algorithms.  The  mirrors  were  driven  by  the  interface 
electronics  package  (the  rack).  In  order  to  drive  the  mirrors,  this 
device  required  digital  information  that  defines  the  frequency  and 
the  voltage  levels  of  the  signals  it  sends  the  mirrors,  and  a 
representation  of  the  nutation  pattern  in  a  form  that  can  be  used  to 
drive  each  mirror. 

The  programs  written  to  control  this  hardware  were  duplicated 
in  a  form  that  could  be  used  by  an  Apple  computer  equipped  with  a 
digital  to  analog  converter.  (These  programs  were  written  in  "Basic" 
and  in  assembly  language  for  the  6502  microprocessor. )  The  Nova 
computer  was  thereby  freed  from  this  chore. 

The  AIT  rack  also  contains  switched  integrators  and  sample 
and  hold  registers  used  in  the  taking  of  data,  These  hardware 
generated  timing  signals  were  used  by  the  Nova  software  to  determine  when 
to  read  the  sample  and  hold  registers. 

The  program  tXPDSK  contains  the  code  necessary  to  take  data 
from  the  sample  and  holds,  scale  the  data,  and  store  the  data  to  a 
disk  file.  EXPDSK  incorporates  the  same  header  and  file  name 
creation  subroutines  as  SIMDSK.  The  files  it  generates  have 
exactly  the  same  format  at  the  SIMDSK  files,  except  for  a  difference 
in  the  letter  code  in  the  file  name. 
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In  order  t.  simplify  data-taking  from  the  experimental  apparatus, 
two  programs  -  variants  of  EXPDSK  and  Image,  respectively  named 
SNAP  and  FSIMAGE  -  were  created.  These  ^wo  programs  are  linked  to 
one  another  by  a  Fortran  "chain",  so  that  each  program  invokes  the 
other  in  a  continuous  process  of  taking  data  and  viewing  the  results. 
Minor  changes  of  the  program  flow  have  been  written  into  these 
programs  to  make  them  more  convenient  to  use  in  this  application. 

There  were  stringent  limitations  imposed  on  the  rate  at  which 
data  could  be  taken.  These  limitations  derived  from  those  of  the 
Data  General  digital  to  analog  converter  hardware.  Scanning  four 
channels,  with  external  sample  and  hold  registers  provided  by  the 
AIT  rack,  no  more  than  3500  sample  points  per  second  could  be 
accommodated.  These  produced  an  uppe*-  limit  of  approximately  5 
frames  per  second  on  the  imaging  rate. 
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8.3  1  he  jGeni sco  Color  D i  splay 

Another  major  software  development  project  undertaken  during 
AIT  Phase  II  relates  to  the  use  of  the  Genisco  Color  Processing 
System.  This  color  processor  permits  a  real  time  display  of  the 
8x8  pixel  field  calculated  by  the  image  processing  equipment. 

Use  of  this  color  processor  required  the  development  of  a  Fortran- 
callable  assembly  language  interface  (PGPDR,  PGPDM)  to  the 
processor  and  also  an  assembly  language  program  to  load  the  vendor 
supplied  color  processor  operating  system  to  the  programmable  color 
processor  (PGP).  Using  the  programming  language  provided  by  the  PGP 
operating  system  and  using  the  Fortran-callable  interface,  code  was 
written  to  create  pseudo-color  displays  of  the  8x8  pixel  image, 
and  other  code  was  written  to  permit  x-y  graphical  displays  to 
appear  on  the  color  monitor.  An  x-y  plotting  utility  permits  a 
6-color  display  of  the  4  detector  vector  waveforms  and  a  measurement 
grid.  Displays  of  other  functions  occur  in  other  programs.  The 
fast  and  easy-to-use  graphing  and  imaging  tools  developed  for 
this  contract  will  be  used  again  and  again  in  future  work. 

With  regards  to  the  real  time  display  of  the  8  x  8  pixel  image, 
the  Genisco  system  was  originally  chosen  for  its  fast  data  channel 
interface  capability  with  Data  General  equipment,  and  its  ability 
to  produce  image  frames  at  high  rates.  A  frame  rate  of  approximately 
16  frames  per  second  was  accomplished,  using  the  relatively  high 
level  PGP  operating  system  executive  language  provided  by  Genisco. 
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Using  a  lower  level  Geriisco  color  graphics  processing  assembly 
language,  further  improvements  in  frame  rates  can  likely  be 
accompl ished. 

These  imaging  and  graphing  subroutines  (using  the  Genisco 
PGP  system)  have  replaced  the  terminal  and  printer  graphing  routines 
previously  used  by  Image  and  SIMDSK.  These  subroutines  and  their 
variants  have  also  been  used  extensively  in  other  programs. 
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8 •  4  General  Disk  File  Handling  JJ t.iliti e  s 

A  general  disk  file  handling  package  developed  for  the  AIAO 
contract  was  adapted  for  use  in  AIT  software  and  was  widely  used. 

The  package  includes  routines  to  write  out  to  disk  a  parameter 
header  and  sets  of  data  and  other  routines  to  read  back  the  header 
and  the  data  sets.  There  is  also  a  file  naming  utility  and  an  error 
handling  utility.  The  naming  utility  creates  a  disk  file  with  a 
unique  name  based  on  the  source  of  the  data  (optical  bench,  or 
simulation,  or  other)  and  also  derived  from  the  date  and  time  of  the 
file  creation.  The  error  handler  sends  error  messages  to  the 
console  and  returns  control  of  the  program  to  ari  appropriate  point 
in  the  program  (that  is,  it  provides  an  error  return  as  well  as  a 
normal  return). 

These  routines  were  designed  with  the  goal  of  being  rriude 
general  enough  to  accommodate  data  sets  of  different  sizes  and 
structure,  for  example,  with  different  nutation  patterns  and 
sampling  schemes.  But,  they  were  made  modular  enough  and  simple 
enough  to  allow  easy  implementation  and  adequate  standardization 
for  uniformity  of  use. 
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9 . 0  EXPERIMENTAL  BREADBOARD  RESULTS 

The  experimental  breadboard  system  was  set  up  for  visible 
wavelength  operation  and  interfaced  to  the  DG  Nova  through  the  All 
Interface  Electronics  rack  and  the  DG/DAC  laboratory  interface.  A 
binary  intensity  bar  target  mask  was  installed  n  the  target 
projector,  and  the  nutation  amplitude  adjusted  so  that  the  bar- 
presented  as  having  dimensions  of  approximately  2x4  pixels  in  the 
reconstruction  field.  Because  of  the  anamorphic  field  lens, 
the  pattern  scanned  had  a  roughly  5:3  aspect  ratio  resulting  in  the 
distorted  reconstruction  field  of  Figure  9.1.  As  a  result,  the  bar 
would  appear  to  be  about  2x4  pixels  in  extent  for  one  orientation 
but  closer  to  1  x  6  pixels  when  rotated  90°. 

If  required,  the  aspect  ratio  can  be  corrected  either  by  an 
anamorphic  re-imaging  system  or  a  simple  deviation  prism-corrector 

The  nutation  drive  was  set  for  a  frame  rate  of  .5  Hz.  This 
slow  speed  was  found  desirable  to  minimize  the  phase  shift  in  the 
linear  galvo  response.  Iri  this  experiment,  no  direct  sensing  of 
the  galvo  mirror  positions  was  used,  ei titer  for  mirror  control  or 
system  synchronization.  In  the  high  speed  system,  mirror  position 
is  sensed  dynamically  and  used  in  a  feedback  loop  to  stabilize  the 
nutation  pattern;  here  the  phase  shift  pattern  does  not  arise. 

Using  the  SNAP  program,  data  sets  were  taken  from  the  bench, 
signal  conditioned,  and  converted  to  output  images.  Systematic 
phase  offsets  were  removed  using  the  shifting  algorithms  and 
DC  restoration  was  performed  on  the  experimental  waveforms,  as 
described  in  Section  4. 
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A  comparison  between  the  simulated  detector  waveform  for  a 
horizontal  bar  and  the  experimental  waveform,  as  measured  at  the 
output  of  the  switched  integrator,  is  shown  in  Figure  9.2.  The 
lower  trace  is  the  simulation  of  quadrant  1,  while  the  upper  is 
quadrant  1  experimental. 

The  results  of  reconstruction  of  the  experimental  waveform  are 
shown  iri  Figure  9.3.  Figure  9.3a  is  the  experimental  result,  and 
9.3b  is  the  simulated  reconstruction,  reproduced  for  comparison. 

The  experimental  bar  is  imperfect  but.  good;  much  of  the 
reduction  in  uniformity  from  simulation  is  due  to  the  slightly 
smaller  size  of  the  real  oar  and  the  fact  that  the  bar  is  not 
precisely  registered  on  the  reconstruction  field  The  rendition 
i c  accurate,  and  is  the  final  proof  of  the  efficacy  of  the  entire 
AIT  imaging  approach. 

The  images  in  Figure  9.4  a.  -  h  are  taken  from  a  sequence  of 
32  reconstructions  of  the  bar  as  it  was  rotated  through  360°.  The 
display  primarily  shows  the  different  effects  of  coarse  image 
sampling.  In  Figure  9.4  ,  the  end  of  the  bar  which  is  rotating 
eccentrically  actually  touches  the  edge  of  the  reconstruction  field; 
the  results  of  aliasing  of  outside  power  are  not  apparent,  indica¬ 
ting  that  the  algorithm  is  robust  with  respect  to  small  nurtures 
or  the  image  from  the  recons truction  definition  field. 

The  reconstruction  o 17  images  from  real  input  data  using  the 
AIT  breadboard  optics  was  the  major  goal  of  the  AIT  program.  It  is 
now  possible  to  carry  out  a  series  of  thorough  experimental  studies 
using  the  AIT  system. 
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li 0 ,.  0  SUMMARY  AND  CONCLUSIONS 

The  AIT  program  has  in  total  been  highly  successful.  The 
basic  goal  of  the  program  has  been  reached:  namely,  a  new  form  of 
imager-tracker  in  which  high-accuracy,  high-update  rate  centroid 
tracking  is  obtained  simultaneously  with  medium-speed  imaging  has 
been  demonstrated.  In  the  process  a  number  of  significant 
theoretical  and  technological  contributions  have  been  made: 

1)  The  AIT  class  of  imaging  algorithms  has  been  analyzed 
and  formally  characterized.  A  general  method  for  obtaining  image 
reconstruction  algorithms  for  any  resolution  has  been  defined. 
Requirements  for  spatial  sampling,  computational  accuracy,  and 
signal  conditioning  have  been  established.  The  accuracy  of  image 
reconstruction  has  been  quantitatively  examined  in  simulation.  The 
signal -to-noise  characteristics  of  the  imaging  process  have  been 
defined,  and  the  results  compared  with  other  imaging  techniques. 
Particular  image-enhancement  characteristics  of  the  reconstruction 
process  have  been  discovered. 

2)  A  class  of  tracker  algorithms  suitable  for  the  AIT 
process  has  been  defined.  A  prescription  for  generating  an  algorithm 
appropriate  to  any  object  size  at  any  resolution  has  been  obtained. 
The  tracker  algorithm  has  beer!  extensively  characterized  in  simu¬ 
lation  and  found  to  have  excellent  properties,  displaying  good 
performance  far  outside  its  optimal  dynamic  range.  The  tracker 

S/M  performance  has  been  estimated  and  approaches  that  of  an  ideal 
quad-cel  1  tracker. 
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3)  A  powerful  new  digital  signal  processor,  the  PMP,  has 
been  designed,  constructed,  and  has  successfully  executed  both 
imaging  and  tracking  algorithms.  Tracking  computation  speed  was 
demonstrated  at  a  full  8  kHz  rate,  while  image  rate  reached  a 

speed  of  2  seconds/frame.  A  fast  vector  processing  unit,  the  VMM, 
was  designed  and  constructed  to  exceed  the  30  Hz  frame  rate  goal, 
but  was  not  tested. 

4)  A  portable  AIT  optical  breadboard  was  designed  and 
fabricated,  and  used  to  demonstrate  reconstruction  of  real  images. 
The  system  is  capable  of  visible  and  infrared  operation,  and 
includes  a  dynamic  target  simulator  with  control  of  target  rotation, 
position,  and  contrast.  The  breadboard  can  operate  at  variable 

low  nutation  rates  with  the  DG  Nova  computer  or  at  fixed  high  rate 
with  the  PMP.  A  special  4-inirror  scanner  technique  was  invented 
and  its  properties  formally  characterized,  so  that  an  appropriate 
nutation  pattern  can  be  defined  for  any  resolution. 

5)  A  large  body  of  experimental  software  was  written, 
including  codes  for  simulation  of  tracking  and  imaging,  experimental 
data  collection,  PMP  development,  and  color  displav  support.  The 
software  is  in  the  form  of  many  utility  routines  called  by  a  few 
large  programs,  allowing  flexibility  in  configuring  the  programs 
for  new  tasks.  The  AIT  code  development  forms  a  substantial  basis 
for  conducting  new  investigations  into  rel.ted  forms  of  imaging 

and  signal  p ro c. es s  i ng . 
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With  the  AIT  breadboard  system  complete,  experimental 
characterization  of  the  imaging  and  tracking  algorithms  may  be 
carried  out.  This  would  be  of  particular  interest  in  the  infrared 
mode  of  operation.  Other  types  of  detectors  can  be  employed,  since 
the  algorithm  formulations  are  quite  general  and  not  limited  to 
use  with  the  quad  cell  detector  configuration.  The  image  processing 
properties  of  the  technique  can  also  be  explored  since  the  appro¬ 
priate  software  and  hardware  tools  are  in  place. 

The  AIT  program  output  has  met  or  exceeded  most  of  its 
goals  and  expectations,  and  can  likely  form  a  nucleus  for  significant 
f i . .  ther  developments. 
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APPENDIX  I:  SOFTWARE  SUMMARY 

This  appendix  contains  a  compilation  of  all  the  software 
developed  for  the  AIT  program.  The  documentation  is  organized 
approximately  along  the  lines  of  usage  within  tiie  AIT  program. 
Main  programs  are  listed  first,  and  all  subroutines  last.  This 
listing  is  intended  as  both  an  overview  of  the  code  and  as  a 
guidebook  for  new  users  of  the  system.  Further  information  is 
available  in  "help"  files  located  on  the  disks  where  the  source 
files  are  stored,  and  in  the  more  general  support  documentation 
for  the  Data  General  and  Genisco  equipment  used. 
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General  purpose  TTY  plot  subroutine. 


Used  by  PlTUV  (See  PLTWV.FR).  This  routine  plots  the  data  passed  to  it  by  the  program 


giving  dimensions  (rows  first)  and  data  following  (row  index  varying  fastest).  Either 
read  or  double  precision  will  work  {X  switch  on  compile  for  double)  --  but  not  mixed 
files  without  d:rrcnsions  at  the  beginning  will  also  work  as  long  as  the  calling  routine 
positions  channels  at  the  start  of  the  data  area. 
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-PGP  --  Continues  write  of  graphics  file  to  PGP  fron  Nova 
DPGP  --  Scads  oata  from  PGP  into  Nova  memory. 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  playis  and  executes  r esecorch ,  development,  tes t  and 
selected  acquisition  programs  In  support  of  Command,  Co yit.ro l 
C c  .muyu.catio ns  and  Intelligence  {Ch)  activitieA.  Technical 
and  engineering  support  within  areaA  of,  technical  competence, 
it  provided  to  ESV  Program  Offices  ( POa )  and  other  BSD 
elements.  The.  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  o.nd  control,  sur- 
v  ell  lane  e  of  ground  and  aerospace  obje.cts,  intelligence  data 
collection  and  handling,  information  system  technology 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliaoilitij,  maintainability  and 
compatibility. 
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